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Abstract 

This article details the computation of the two-point correlators of the convergence, E- and B- 
modes of the cosmic shear induced by the weak-lensing by large scale structure assuming that the 
background spacetime is spatially homogeneous and anisotropic. After detailing the perturbation 
equations and the general theory of weak-lensing in an anisotropic universe, it develops a weak 
shear approximation scheme in which one can compute analytically the evolution of the Jacobi 
matrix. It allows one to compute the angular power spectrum of the E- and B-modes. In the 
linear regime, the existence of iJ-modes is a direct tracer of a late time anisotropy and their 
angular power spectrum scales as the square of the shear. It is then demonstrated that there must 
also exist off-diagonal correlations between the E-modes, .B-modes and convergence that are linear 
in the geometrical shear and allow one to reconstruct the eigendirections of expansion. These 
spectra can be measured in future large scale surveys, such as Euclid and SKA, and offer a new 
tool to test the isotropy of the expansion of the universe at low redshift. 
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I. INTRODUCTION 


A. Motivations 

The standard model of cosmology describes our universe with a very simple solution 
of general relativity describing a spatially homogeneous and isotropic spacetime, known as 
the Friedmann-Lemaitre solution. It is assumed to describe the geometry of our Universe 
smoothed on large scales. Besides, the use of the perturbation theory allows one to under¬ 
stand the properties of the large scale structure, as well as its growth from initial conditions 
set by inflation and constrained by the observations of the cosmic microwave background 
(CMB). It is a very successful model and allows one to deal with all existing observations 
in a consistent way with only 6 free parameters [Tj from primordial nucleosynthesis (BBN) 
to today, involving mostly general relativity, electromagnetism and nuclear physics, that 
is physics below 100 MeV and well under control experimentally (see, e.g., Refs. PHI] for 
standard textbooks). 

The construction of the cosmological model depends on our knowledge of micro-physics 
but also on a priori hypothesis on the geometry of the spacetime describing our universe. It 
relies on 4 main hypothesis (see Ref. [5] for a detailed description): (HI) a theory of gravity, 
(H2) a description of the matter and the non-gravitational interactions, (H3) symmetry 
hypothesis, and (H4) an hypothesis on the global structure, i.e. the topology, of the universe. 
The hypothesis HI and H2, that refer to the physical theories, are not sufficient to solve 
the field equations and we need an assumption on the symmetries (H3) of the solutions 
describing our universe on large scales. 

Among the generic conclusions of this standard model is the need of a dark sector, includ¬ 
ing dark matter and dark energy, which emphasizes the need for extra degrees of freedom, 
either physical (new fundamental fields or interactions) or geometrical (e.g., a cosmological 
solution with lower symmetry). This has driven a lot of activity to test the hypotheses of the 
cosmological model. In that debate, weak-lensing is a key observation to test general rela¬ 
tivity on cosmological scales jH] and to constrain the scale on which the fluid limit holds [7]. 
It complements tests of the other hypothesis such as the equivalence principle [5] and the 
Copernican principle [9]. Our first motivation is thus to provide a new test on the isotropy 
of the expansion at late time, hence providing a new test of the standard cosmological as¬ 
sumption. Any detection of a violation of a symmetry of the background spacetime would 
have important implications in terms of model building and on the understanding of the 
dark sector. 

While in the standard ACDM model the cosmological constant A is the source of the 
acceleration of the universe, many models have been proposed to explain the acceleration 
of the cosmic expansion. The property of the dark sector is often modelled as a fluid 
with an equation of state, Pde = ^Pde, relating its pressure to its energy density. Such a 
phenomenological parameterization allows to characterize the ability of different surveys 
to actually demonstrate that w = —1, as expected for a cosmological constant. Among 
the plethora of dark energy models, many enjoy an anisotropic pressure H) and thus may 
trigger a phase of anisotropic expansion at late time when dark energy starts influencing the 
dynamics of the universe. This is for instance the case of magnetised dark energy [TO], ITT] , 
solid dark matter mm induced by a network of frustrated topological defects, bi-gravity 
models [QE], anisotropic dark energy mm and in models in which the backreaction mi 
of the large scale structure on the background evolution is the source of the acceleration. 
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This has led to the development of a phenomenological parameterisation of the anisotropic 
pressure in terms of an anisotropic equation of state as II* = Ain*p de [T8H2T] . Our second 
motivation is thus to propose new observational tests on the anisotropic pressure of the 
dark energy sector, hence constraining another phenomenological deviation from a pure 
cosmological constant. 

When concerned by anisotropic expansion, we can distinguish between 2 classes of models, 
that allows one to divide the different methods to constrain anisotropy. Remind that any 
perturbed quantity, A" say, such as the gravitational potential, the density contrast, etc. 
can be split, in Fourier space, as the product of an initial configuration and a transfer 
function as X(t,k) = T\{t, fc)Xj(fc). First, early anisotropic models (such as anisotropic 
inflation) have anisotropic initial conditions (in the sense that the correlation functions 
of the initial perturbed quantities is such that (Xi(k)X*(k')) ^ Px{k)5(k — k ')) while 
the transfer functions are independent of direction (i.e. Tx(t,k) = T\(t,k )) because the 
geometry has isotropized at later times. Second, late time anisotropic models have been 
isotropic during most of the history of the universe (hence enjoying isotropic correlation 
functions, e.g. (Xi(k)X*(k')) = Px(k)8(k — k')) while their transfer function at late time 
is anisotropic, i.e. Tx(t, k) ^ T\(t, k). These two types of models have a huge difference in 
the way one attacks observational constraints. In particular the propagation of light is only 
affected in the second class of models. 

Without any source during inflation, any primordial anisotropy is washed out [221 22] 
by the expansion. It was however demonstrated that it affects the construction of the 
Bunch-Davies state ra so that it lets very specific signatures on the primordial power 
spectrum [24] 126] and affects the onset of inflation 123 . Such deviation from isotropy can 
be constrained by CMB observations [[281435] . An early, post-inflationary, anisotropy also 
affects the synthesis of light elements during primordial nucleosynthesis [36] (mostly because 
it affects the expansion rate). 

Tests of a late time anisotropy have mostly focused on the Hubble diagram from type la 
supernovae [21, I57M^5] . An anisotropic expansion will influence the transfer function so 
that it can also be constrained by the study of the large scale structure [4TH55] and of the 
CMB [561 loll] . It was argued that supernovae data leads to Aw < 2.1 x 10 -4 [H], and that 
next generation galaxy surveys are capable of constraining anisotropies at the 5% level [21j 
in terms of the anisotropic equation of state. 

In this article, we follow our former analysis [60] on the imprint of a late time anisotropy 
on weak-lensing. According to the standard lore pH], in a homogeneous and isotropic 
background spacetime, weak-lensing by the large scale structure of the universe induces a 
shear field which, to leading order, only contains R-modes. It was demonstrated in Ref. [60] 
that, even in the linear regime, anisotropic expansion will reflect itself in the existence of 
non-vanishing R-modes. The level of R-modes is used as an important sanity check during 
the data processing. On small scales, R-modes arise from non-linear effects [62] intrinsic 
alignments [63] . Born correction, lens-lens coupling [63], and gravitational lensing due to 
the redshift clustering of source galaxies [65] . On large angular scales in which the linear 
regime holds, it was demonstrated |60j that non-vanishing R-modes would be a signature 
of a deviation from the isotropy of the expansion, these modes being generated by the 
coupling of the background Weyl tensor to the R-modes. 

While light propagation in strictly homogeneous Bianchi universes has been widely in- 
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vestigated |66j |67], the analytic computation of the Jacobi matrix was only determined re¬ 
cently [68] (see also Ref. [69]). This article focuses on the computation of the Jacobi matrix 
taking into account cosmological perturbations at linear order in a spatially homogeneous 
anisotropic Euclidean spacetime of the Bianchi / family. We provide all the technical tools 
(perturbation theory, light propagation, expression of the observables). The application of 
our formalism is exposed in our companion paper 123] in which we compute the expected 
signals for the Euclid [71] and SKA [72] observations. 

Among our main results, we emphasize that, as soon as local isotropy does not hold 
at the background level, there exist a series of weak-lcnsing observables that allow one to 
fully reconstruct the background shear and thus test spatial isotropy. More precisely, as a 
consequence of the non-vanishing of the E-modes, it can be demonstrated that 

1. the angular correlation function of the E-modes, C BB , is non vanishing [611], and scales 
as the square of the ratio of the geometric shear to the Hubble expansion rate, a 2 /'H 2 \ 

2. the E-modes correlate with both the E-modes and the convergence k leading to the 

off-diagonal cross-correlations {Bi m E^ ±lm _ M ) and in which E £m and 

Bi m are the components of the decomposition of the E- and E-modes of the cosmic 
shear in (spin-2) spherical harmonics and the components of the decomposition 
of the convergence in spherical harmonics. These two correlators scale as 

3. the deviation from isotropy also generates off-diagonals correlations among k and E 

modes, (E £m E^ ±2m _ M ), and (E £m K* i±2m _ M ). These three correlators 

scale as <7/7/; 

4. for each type of correlator, there are five values of M so that in principle they can be 
used to reconstruct the five components of the geometric shear aij. 

This last point is very important since it exhibits a rigidity between independent observables 
that can be used to control systematic effects. 


B. Structure of the article 


Section [TT] summarizes the description of the spacetime at the background level (§ II A) 


and for linear perturbations (§ IIB). For the sake of clarity, the theory of gauge invari¬ 


ant perturbations is detailed in Appendix [A} It also introduces the parameterisation of an 
anisotropic dark energy sector. The main variables required to describe the evolution of the 
background spacetime are summarized in Table | 


Section III describes the propagation of a light bundle (§ III A) and presents in § III B 


the central equation for our analysis, namely the Sachs equation 


d 2 

_ T> a — 7?“T) C 

dv* b ' 


b j 


for the 2x2 Jacobi matrix V ab , the decomposition of which is presented in § E C] I t 
concludes by specifying these general results to the case of a Bianchi I spacetime (§ HID), 
focusing on the technical but useful use of a conformal transformation. The main variables 


required to describe the evolution of a geodesic bundle are summarized in Table IB 


Since the geometric shear is obviously small, we develop in Section [TV] an approximation 
scheme referred to as small shear limit in which a/7~L is considered as a small parameter. 
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Symbol 


Meaning 


Appears at Eq.: 


Ab ”, ■ ■ ■ 
i,3, ■ ■ ■ 
hh ■ ■ ■ 
0 


a 

H 

Pi 


a ij 

a ij 

#0 

®(f 

Aw J - 

Pa 


Formal space-time indices 

Cartesian spatial indices 

Spatial tetrad indices 

Time tetrad index 

Average scale factor 

Cosmic time Hubble expansion rate 

Log of directional scale factors 

Geometrical (cosmic time) shear 

Geometrical (conformal time) shear 

Timelike vector of background tetrad 

Spacelike vector of background tetrad 

Timelike vector of perturbed tetrad 

Spacelike vector of perturbed tetrad 

Equation of state of dark energy anisotropic stress 

Homogeneous perturbation of the Euclidean metric. 

TABLE I: Table of most used quantities describing the 
background spacetime. 


( 2 . 11 ) 

(2.13) 

M 

|2T8a| 

(P|) 

(2A) 

2J3) 

( pi 

( 2 . 11 ) 

MM 

T6) 

(p6|) 

S) 


Symbol 

Meaning Appears at Eq.: 


Null geodesic tangent vector 

(3.1 


& 

Conformally null geodesic tangent vector ( 

3.33 

) 

z 

Redshift 

(3.4 


n Q 

Initial observed direction 

(3.5 


n a 

Sachs basis 

(3.8 


n± 

Helicity basis ( 

3.11 

) 

Va 

Components of the connecting vector in the Sachs basis ( 

3.13 

) 

7t a b 

Optical tidal matrix ( 

3.13 

) 

Dab 

Jacobi matrix ( 

3.14 

) 

Dab 

Conformal Jacobi matrix ( 

3.37 

) 

Da 

Background angular diameter distance ( 

3.17 

) 

Da 

Angular diameter distance ( 

3.19 

) 

K 

Convergence ( 

3.17 

) 

'Jab 

Cosmic shear ( 

3.17 

) 

V 

Rotation ( 

3.17 

) 


Multipolar coefficients of the convergence ( 

3.23 

) 

Vim 

Multipolar coefficients of the rotation ( 

3.23 

) 

Dim 

Multipolar coefficients of the cosmic shear E-modes ( 

3.25 

) 

Dim 

Multipolar coefficients of the cosmic shear E-modes ( 

3.25 

) 


TABLE II: Table of most used quantities describing the 


propagation of a geodesic bundle. 







































































































{ 0 , 1 } 


x« 0 ’°>n° n° 

x< 0 ’ 0 > n° 



>.°> 



FIG. 1: Comparison of the geodesic in the approximation at order {n,p}. In order to adopt an 
observer based point of view, we need to relate the local direction of propagation n^ n,p ^ to n°. 
The transports for the 3 orders of perturbations are respectively detailed in Eqs. ( 6.12[ ) and (6.37) 
that determine xd n,p ^ (x, n°) that can be further split in a radial component Sr^ n,p ^ (y, n°) and an 
orthoradial contribution that defines the deflection angle a a ^ n,p ^ {xi n °)■ 


We then use a two-parameter expansion scheme in which both a/7-L and the perturbations 
of the metric, say <f>, are small. Thus, a given order {n,p} corresponds to term of order 
(cr/"H) n $ p . In this approximation, the structure of our computation is the following. We 
start from the fact that the Sachs equation can be rewritten as (see Eq. (4.34)) 


d 2 V ab + 1 dk°dV ab 


^IZacT), 


ac^cb i 


d \ 2 ' k° dx dx ( k °) 2 

where x is t he co ordinate along the lightcone in the background Friedmann-Lemaitre space- 


time (see 
the form 


(see § IVD). At order {0,0}, 7£^’°* = 0 and fc 0 ! 0 ’ 0 ! = —1 so that the equation takes 


d 2 D { ”’ 0} 

-^— = 0 

d X 2 


and can be integrated trivially (see § VIB). We then expand this equation order by order 
so that it formally takes the form (since fc| 00 } = — 1) 



dx 2 


g{n,p} 


in which the source term contains contribution from lZ ab and k° up to order {n,p}, and from 
T> ab at lower order. The effects to be taken into account are then 

1. the tensor and vector contributions to 7 Z ab: which starts at order {1,1} and the con¬ 
tribution of the scalar modes at the relevant order; 

2. the evolution of all the perturbative modes, that is of the transfer functions, which is 
decomposed as 

T x (k,t)=Tf»(k,t)+Tp 1 \k t t) 

since the order {0,0} and {1,0} correspond to homogeneous solutions. This requires 
to solve the equations of Appendix |A| 

3. In order to determine k°, we also need to solve perturbatively the geodesic equation. 
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4. A source observed in direction n° at distance x, is located at a spacetime point P{ n , p } 
and its contribution depends on the local direction of the tangent vector to the geodesic 
in which determines the local Sachs basis in P{ n , p }- We shall thus proceed with 

2 operations: 

(a) transport P {ritP} to P {0 , 0 } (see Fig. [l|, 

(b) transport ^{ n , P } and the local Sachs basis (see Fig. [2]). 

This is what we call the central geodesic approximation and the possibility to go beyond 
this approximation is sketched in Appendix [Cj We however stick to this approximation, 
which is sufficient in the small shear approximation. At lowest order, it corresponds 
to the usual Born approximation but at higher order there are post-Born corrections 
to include. 


Physical geodesic Fiducial FL reference geodesic 



n< 0 ’°> 


n° 

< 

►-♦ 


{0,0} 



FIG. 2: The local Sachs basis at a point P{ n , p } on the geodesic has to be transported to the 
point P {0 ,o} on the fiducial null geodesic of the background Friedmann-Lemaitre spacetime. This 
implies to perform a transport on the tangent space at the same time that the point of observation 
is transported in real space. The transports for the 3 orders of perturbations are respectively 


detailed in Eqs. ( |6.15[ ), ( |6.34[ ) and ([6434]) for n and in Eqs. ( |6.17[ ), ( |6.35[ ) and ( |6.54[ ) for the Sachs 
basis. 


Section [V] describes the computation of the angular correlation. Our philosophy is to 
adopt an observer point of view, that is, to compute all quantities on the celestial sphere of 
the observer. Given the previous perturbative expansion scheme, any observable X s of spin 


s can then be formally expressed as (see Eq. (5.4)) 




X s {x 3 , n o)m s 0 = / S S (x s ,x,n 0 )m s 0 dx, 


for a source S x ‘ located at y s and observed in direction n°. According to the spin s of 
the quantity we can expand in the proper spherical harmonics with respect to n°. This 
allows one to define the expansion of all the quantities in term of spherical harmonics. As 
a byproduct, we demonstrate in § VC that the 5 off-diagonal correlators {B £rn E^ ±lm _ M ) , 


(B im K* l±lm _ M ) , (E £m E* ±2m _ M ), and {E £m K* l±2m _ M ) are non-vanishing. 

Equiped with all these tools, we compute these correlators and the angular power spec¬ 
tra of the E- and 5-modes in Section VI order by order. Order {1,0} recovers the non- 
perturbative analysis of Ref. [68] while order {0,1} recovers the standard case of lensing 


10 















Symbol 


Meaning 


E 

S a 

S a b 

D a 


im 


^2m 

s± 

m o 

S Xs 

rr-iJCg 

± lm 
Xg-rLM 
± tm 

B x 

Cf E 

cf B 

B 

a a 

w a 

V 

M(x) 

m 


Scalar shear on the sphere 
Vector shear on the sphere 
Tensor shear on the sphere 
Covariant derivative on the sphere 
spin-raising and lowering operators 
Multipole of the scalar shear 

First derivative of the scalar shear in the helicity basis 
spin-s polarization basis 
Source of the held X s 

Anisotropic transfer function of the held X s 
Multipoles of the anisotropic transfer function 
Multipoles of the E- and B-modes of the held X s 
Angular power spectrum of the B-modes 
Angular power spectrum of the B-modes 
Scalar perturbation of the spatial metric 
Deflection angle 

Perturbation of propagating direction 
Delhection potential 
Source distribution 
Primordial power spectrum 


Appears at Eq.: 

( 4.16 ) 

( 4.16 ) 

( 4.17 ) 

( 4.21 ) 
( 4.25 ) 

(4.27) 

(pf 

5 A) 

( pi 

( 5.15 ) 

( 5.17 ) 

( 5.21 ) 
(5.21) 

s 

( pi 

( 6.17 ) 

( 6.27 ) 

(£43) 

(5.14) 


TABLE III: Table of most used quantities describing the 
propagation of a geodesic bundle. 


by large scale structure in the linear regime; 


VI E| gives all the details of the computa- 

In 


tion at order {1,1}. This allows us to discuss the dominant contribution in Section VII 
particular, we argue that the dominant term for the cosmic shear is given by 


Tab Q? D c D^ a D)^(p 


where a c is the deflection angle at order {1,0}, Lp the deflection potential, and D a the 
covariant derivative on the celestial sphere. 

Many technicalities are gathered in the appendices: linear perturbation theory (§ 0. 
the expressions of the geometric quantities at first order in perturbation that are need to 
compute the source term of the Sachs equation ( § [B| , details on the lensing method (§ 0 
and a catalog of useful mathematical identities (§ |D[). Throughout this work we adopt units 
in which c — 1. 
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II. SPACETIME STRUCTURES 


A. Background spacetime 

1. Geometry 


At the background level, the universe is described by a spatially Euclidean, homogeneous, 
and locally anisotropic solution of the Einstein equation filled with a perfect fluid. Its metric 
takes the general form (see Refs. (73H77] for general references on Bianchi spacetimes) 


3 


i \2 


ds" = -dt z + 2^Xt(t)(dx\ 

i =1 

= —dt 2 + a 2 {t) r ) i j[t)dx' l dx\ 


(2.1a) 

(2.1b) 


where a(t) = yjX\ (t)X 2 (t)X 3 (t) is the volume averaged scale factor and t the cosmic time. 
We define the tangent vector to the fundamental comoving observer by u^dx^ = dt. It 
is normalized such that = —1. The spatial metric y^- and its inverse y*- 7 can be 

decomposed as 

7 *j(t) = exp[2 PiitySij , 7 13 (t) = exp[-2:ft(t)]<5* J , (2.2) 

with the constraint 


X> = o 


(2.3) 


i —1 


that ensures that the comoving volume remains constant (i.e, 7 = Y' ]A /ij — 0). Note that, as 
a consequence of Eq. (2.2), some spatial directions should contract while others grow BS- 
Note also that that there is no sum on i in the definition of 7 \j and Latin indices {i,j, k ... } 
are raised with y *- 7 and lowered with 7 
The geometrical shear is dehned as 


a ij — 2 Tij 


(2.4) 


where a dot refers to a derivative with respect to cosmic time. We shall also use the conformal 
time 7 dehned by a(rj)drj = d t, and denote derivatives with respect to it by a prime. Thus, 
the conformal shear is dehned as 


a ij — 2 Aj a<T ij ■ 


The amplitude of the shear is dehned by 


.2 _ 


cr = a= 


5W 2 

i= 1 


and 


2 _ 


a~ = aija 11 = a 


3 

- 2 E 




l— 1 


2—1 


(2.5) 


( 2 . 6 ) 


2. Decomposition of the geometric shear 

The shear, being a symmetric and traceless spatial tensor (a\ = 0), has 5 degrees of 
freedom, three of which correspond to the Euler angles necessary to express the shear in a 
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general basis. By choosing the Cartesian basis (2.2), we have set these three angles to zero so 
that we are left with only 2 degrees of freedom, namely the three fa with the constraint (2.3). 
The components of the shear can thus be expressed as 


crij{t) = A exp [2fa(t)]S i:j , a 13 {t) = fa exp[-2fa(t)]5 lJ , a)(t) = fad) . (2.7) 

These two independent degrees of freedom can also be decomposed as a magnitude and an 
angle ip. The first is related to the scalar shear while the choice of the angle defines which of 
the spatial directions are initially expanding. These variables can be obtained by performing 
the decomposition 

fa{t) = awfa ), ( 2 . 8 ) 

with the 3 constants Ci given by 

Q = Sin (</? + Y«) , ie {1,2,3} (2.9) 

and where S is constant. This parameterization obviously satisfies the required constraints, 

^a = o, Yl c ? = s *- ( 2 - 10 ) 

i i 

Therefore, once the Cartesian basis is chosen, we can choose the two constants (<p,S) to 
describe the two degrees of freedom of the shear since a 2 = (SW) 2 . 


3. Spatial triad 

It is convenient to introduce a spatial triad - a set of three orthonormal vectors and co¬ 
vectors; the normalization being defined from 7 \j and fa 3 - related directly to the Cartesian 
coordinates x l . Their components in the coordinates basis read 

= exp [-fa(t)]5 3 i , fa-j = exp [fa(t)]5] . (2.11) 

In such a triad basis, the shear components take the simple form 

&ij = a 1 - 3 - = fadij Oij = a 1 - 3 - = faSij. ( 2 . 12 ) 

Thus, this triad can easily be extended to a tetrad by using the observer’s 4-velocity as the 
normalized time-like vector 


K = % = < = ^ = (2.13) 

4■ Description of matter and field equations 

Concerning the matter sector, we assume it is composed of a pressureless matter fluid 
and a dark energy component. The dark sector is then described by a fluid whose energy- 
momentum tensor enjoys a non-vanishing anisotropic stress, 

Tl = (p + P)u' i u v + P5Z + IP„. (2.14) 


13 




The anisotropic stress tensor is symmetric (n M „ = n„ M ), traceless (11^ = 0) and transverse 
(u^ = 0) which means that it has only 5 degress of freedom encoded in its spatial part 
II ij. Unless we define a microscopic model, we need to use an equation of state for II*. We 
decompose it as 

n; = p de Aw) , (2.15) 

so that the pressure tensor takes the general form 

Pi = pde (w8j + Awl) , (2.16) 

where w is the equation of state relating the isotropic pressure to the energy density and 
Aw\ an equation of state for the anisotropic pressure. From a phenomenological point of 
view, this corresponds to an extension of the dark energy sector, similar to the ansatz (1) 
of Ref. [TT|, which will allow us to address the question of the deviation from the standard 


cosmological constant reference {i.e. w — — 1 and Awl = 0)- 
Defining the Hubble expansion rate by 

H = a/a, 

(2.17) 

the background equations |22j take the form 

3 H~ — K(p m + Pde) + 2^ 2 J 

(2.18a) 

(<j))' = —3 Ha) + kU) . 

(2.18b) 

Pm = 3 H p m , 

(2.18c) 

Pde = —3H(1 + w)p de - UijW 3 . 

(2.18d) 


The first equation is the analogous of the Friedmann equation in presence of a spatial shear, 
the second is obtained from the traceless and transverse part of the Einstein equation and 
dictates the evolution of the shear. The last two equations are the continuity equations for 
the dark matter (P = IT) = 0) and dark energy sector. We have set 1 k = 8ttG = Mp 2 . 


5. Dynamics 


The set of equations (2.18) can be formally integrated. As usual, the dark matter energy 
density scales as 

' a °^ 3 (2.19) 


Pm ~ P mO \ — J 


Equation (2.18b) has a first integral given by 


cr- 


/ aoV 


C* + k / IE I -1 A{a/ao) 


a o 


H 


( 2 . 20 ) 


1 In order to easily check the homogeneity of the equations, we recall that 


[H}~M P , \p\~M%, [k] ~ M ~ 2 , [II] ~ Mp 
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where C* is a constant tensor representing the decaying mode of the shear. Note that if the 
term proportional to C* is not negligible then the shear is not proportional to the anisotropic 
stress so that and fl^ cannot be diagonalized in the same basis. Integrating Eq. (2.18d) 
leads to 

-3(1+id) 


Pde 


Pdeo - / ajU) ( - 


a \ 2+3w d(a/ao) 


H 


( 2 . 21 ) 


or, if one uses the decomposition (2.15), as 

-3(1+™) 


Pde — PdeO 


exp 




( 2 . 22 ) 


In the particular case where w — — 1, this latter equation teaches us that the dark energy 
density does not remain constant. 


B. Linear perturbations 

Cosmological perturbation theory around a Bianchi / background spacetime, in the 
Bardeen formalism, was first investigated in Ref. [ 22l f24] . The perturbed spacetime has 
a metric of the form 


ds 2 = a 2 [— (1 + 2A)dr] 2 + 2Bidx z d'p + ( 7 ^ + hij)dx l dx j ] 


(2.23) 


where A is a free scalar function, Bi = diB + Bi and h VJ = 2 C ( 7 ^ + + 2didjE + 2d(iEj) + 
2 Eij defined together with the usual transversality and trace-free conditions d^B 1 = 0 = 
diE\ E\ = 0 = diE lj . 

As summarized in § A 2 one can define 2 scalar (<h and T), 2 vector (<3>j) and 2 tensor 
(E^) degrees of freedom, which are gauge invariant; see Eq. (A25). Similarly, one can define 


gauge invariant variables for the matter sector, leading to 4 scalar variables (Sp, 5P, v and 
tt s respectively for the density, pressure, velocity and anisotropic stress), 4 vector var iables 


(v l and nY ) and 2 tensor variables (if fA the expressions of which are gathered in Eqs. (A33) 


and (A35). 


Appendix |A| summarizes all the techniques and results needed to study the perturbations, 
including the definition of the Fourier transform (§ | A1), and the construction of the gauge 


invariant variables (§ A 2). It then derives the full set of Einstein equations (§|A3) and the 


conservation equations (§ A4). 


Among the important features that differ from the standard perturbation theory around 
a Friedmann-Lemaitre spacetime, let us mention 

• the fact that only the components /q of the wave-(co)vector are constant so that both 
k l and k are time dependent - see e.g. Eq. (A2); 


the fact that the scalar-vector-tensor modes do not decouple; 

the fact that, even at late time, the two Bardeen potentials are not equal because of 
the anisotropic stress. 
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III. WEAK-LENSING IN A GENERAL SPACETIME 


This section provides the definitions and equations describing the propagation of a 
geodesic bundle (geodesic equation and Sachs equation) in a general spacetime and in the 
particular case of a Bianchi / universe. 


A. Geodesic bundle 

Weak-lensing is concerned with the deformation of an infinitesimal bundle of light rays 
propagating in curved spacetimes. It is thus related to the geodesic deviation equation. 


1. Geodesic equation 

The central quantity in the geodesic equation, obtained as the eikonal limit of Maxwell’s 
equations, is the wave-vector of an electromagnetic wave, k^{y) = dx M /du, where v is an 
affine parameter of a given geodesic x^(v) and defined such that it is zero at the observer 
and increases towards the source. We shall be working in the eikonal approximation where 
kv is a null vector satisfying the geodesic equation 

rV.F = 0. (3.1) 

If we parameterize the bundle of null geodesics by x**(v, s ), where s is a continuous parameter 
labeling each ray of the bundle, then each ray has a wave-vector given by /c M (v, s) = dx^/dv, 
whereas the vector rf 1 = djds gives the infinitesimal separation between two neighbouring 
geodesics of the bundle. The photon wave-vector can always be decomposed in components 
respectively parallel and orthogonal to vJ 1 as 

#* = CT 1 *** = -«** +n" , (3.2) 

where are the components of the local directional vector n, defined such that 

v/rin = 0 , ri^ri 11 = 1. (3.3) 

Once the geodesic equation is solved, any comoving observer with four-velocity normal¬ 
ized such that v/ l Un = — 1, defines the redshift of a source by 

l + Z{V ’ n ° ] "<k^%' (3 ' 4) 

where v is the affine parameter that specifies the position of the source down the lightcone 
and 

n 0 = n(v = 0) (3.5) 

is the direction of observation. The energy of a photon at a given redshift is 

U(v,n 0 ) = U 0 [l + z{v,n 0 )\, U 0 = (k fl u^) 0 . (3.6) 

By definition, the local spacclike vector n is a function of the affine parameter v and of 
the direction of observation observer n Q , that is the spacclike vector pointing along the line 
of sight. 
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2. Geodesic deviation equation 


A (narrow) light beam is a collection of neighbouring light rays. The behaviour of any 
such geodesic, with respect to an arbitrary reference one, is described by the separation (or 
connecting) vector r /A Asuming that all the rays converge at a given event O (the location of 
the observer), 7/^(0) = 0. The evolution of i^(y) along the beam is governed by the geodesic 
deviation equation 


d 2 77 M 

dv 2 

where R^ ua p is the Riemann tensor. 




(3.7) 


3. Sachs basis and screen space 


For any observer whose worldline intersects the light beam at an event different from O, 
the beam has a non-zero extension, since a priori ^ 0. The observer can thus project 
it on a screen to characterize its size and shape. This screen is by essence a 2-dimensional 
spacelike hypersurface, and chosen to be orthogonal to the local line of sight nA Two such 
spatial vectors required to construct a basis for the tangent space, n a with a = {1,2}, are 
defined by the requirement that 

= S ab , = 0 . (3.8) 

With these definitions we can construct a tensor which projects any geometrical quantity 
on the two dimensional surface orthogonal to n 


S'i_n/ — Qfj,i> T MjjUy n^Tiy . 


(3.9) 


Then, with the help of the orthogonality relations (3.8), this two-dimensional screen basis 
can be parallel propagated along null geodesics as [78] 

(3.10) 


= o. 


A basis satisfying the condition (3.8) and propagated according to Eq. (3.10) is called a 


Sachs basis. It is important to note that the basis formed by the vectors n a is defined up 
to an overall rotation around n Q . We can fix this freedom by introducing a spherical basis 
at the observer (i.e., at v = 0) by demanding that {n°,n° 1 n 2 } = {n“,n^, n°}. With this 
choice, the integration of Eq. (3.10) allows us to define a unique three-dimensional basis 
{n r ,n e ,n^}(n°,ii) at each point along the geodesics; see Ref. [60], Furthermore, it will be 
convenient to define a helicity basis as 


n ± = -h (n„ =F in,) = A ( n i =F '"2).. 


(3.11) 


whose components in the n a basis read 


n% = n ± -n a = =F MJ) 


(3.12) 


and are, by construction, constant. 
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B. Sachs equation 


The screen projection of the connecting vector, p a = 72 ^ 77 ^, represents the relative position 
on the screen of the two light spots associated with two rays separated by 77 ^. Similarly, and 
if we set by convention U Q — 1, 9 a = ( dr] a /dv) 0 represents the angular separation of those 
rays, as observed from O. 

The geodesic equation can be recast as |2J an equation for r) a 


as 


d 2 _rja 

du 2 


= KabV b 


Kab = R^apk u k a n^ 


(3.13) 


lZ a b is the screen projected Riemann tensor which can be split into its symmetric traceless 
part lZ{ab) and its trace part 1Z = !Z a bS ab /2. It is also referred to as the optical tidal matrix. 
Furthermore, thanks to the linearity of Eq. (3.13), one can decompose the connection vector 
on the geodesic to its initial derivative as 


This defines the Jacobi map V a b that satisfies the Sachs equation [21179] 

d 2 

dv 2 » b -Kc»b, 

subject to the following initial conditions: 

r\T) a 

iW) = o, -^(o)=K- 


(3.14) 


(3.15) 


(3.16) 


C. Decomposition of the Jacobi matrix and observables 


The Jacobi matrix entering the Sachs equation (3.15) encodes all the information about 
the deformation of a light beam when propagating through a curved spacetime. This 2x2 
matrix can be decomposed in different ways. 

The usual decomposition is described in terms of a convergence n, a rotation V, and a 
shear 7 a b as 

Rabid 0 = Da{v) [(1 + n)I a b + VCab + Jab] ( 3 - 17 ) 


with 

e ab = 2m^ a n+, 7 “ = 0, (3.18) 

and where screen-basis indices a and b are manipulated with I a b = = 5 a b, that is 

with a two-dimensional Euclidian metric. 

A canonical decomposition was introduced in Ref. [68j as 


V a b{v) = D a (v) 


cos sin 7/1 

exp 

-id 

r 2 

— sin 7/1 cos 7/1 

r 2 

r 1 

-V-- 





rotation cosmic shear 


(3.19) 
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According to this decomposition, the real size and shape of the light source is obtained from 
the image by performing the following transformations: (i) an area-preserving shear (r 1; T 2 ), 
(ii) a global rotation 0, (iii) a global scaling. The latter defines the angular distance as 


D a (v) = ^detV ab (v ), 


(3.20) 


which does not assume a ny ba ckground spacetime and perturbative expansion. On the 
other hand, the definition (3.17) introduces the background angular distances D A . Both are 
related by 

D a (v) ~ D a (v)[ 1 + k(v)] . (3.21) 

As for the deformation of the source shape, it is given by the reduced shear 


'D(ab) ^ 'Yab 
Vdet V ab (1 - k) ' 


(3.22) 


Each one of the above observables are dehned on our past lightcone, and, as such, they are 
functions of n Q and v. The convergence and the rotation are scalar functions, and therefore 
can be expanded in terms of scalar spherical harmonics as 

n(n°,v) = y ^Kt m {v)Y em (n 0 ) , (3.23a) 

£,m 

V(n°, v) = J2 V im (v)Y em (n°) . (3.23b) 

£,m 


The cosmic shear, on the other hand, being a spin two quantity, can be expanded in terms 
of the polarization basis as 

7 ab(n°, v) = Y^ 7 A (™°, v)n*n$ . (3.24) 

A=± 

The coefficients ^ can be further expanded in terms of E- and 5-modes on a basis of spin-2 
spherical harmonics as 


7 ± (n°, v) = ^ [E tm {v) ± i B em (v)\ Y^(n °). (3.25) 

£,m 

It should be stressed that we adopt an observer-based point of view. This means that 
all quantities are expressed in terms of (n°,v). In general, n(n°,v) ^ n°, with the obvious 
exception of, e.g., Friedmann-Lemaitre spacetimes and spacetimes with a local spherical 
symmetry for an observer located at the center of symmetry. Therefore, one of the difficulties 
in obtaining cosmological observables as a function of v, or equivalently as a function of the 
redshift z, lies in the determination of these coefficients. 


D. Particular case of a Bianchi I spacetime 

1. Geodesic equation 

The Bianchi / spacetime enjoys 3 Killing vectors, <9*, that allow one to construct 3 con¬ 
served quantities, g(di, k) = ki , along any geodesic. It implies that 

ki = cst, (3.26) 
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on each geodesic so that 


k { = 



k being a null vector, one concludes that co 2 = (k 1 ) 2 = g^kjkj with 


(3.27a) 


co = - f 
a 


£>-***)*. 

\ i = 1 


(3.28) 


It follows that the components of the direction of observation vector are given by 

7ii = ki/co, n l — k l /oo. (3.29) 


The constants of motion ki are then directly related to the direction in which the observer 
in O needs to look to detect the light signal, i.e. the direction of the source n%. The redshift 
of a source is then given by 


1 + z(n 0 ,t s ) 


ujs 


CLq 

a(t s ) 



e -fths )ki 2 
[e-ftho )k i ] 2 


(3.30) 


It is always possible to choose the normalisation such that a Q — 1 and (3i(t 0 ) = 0, but we do 
not make that choice here. 


2. Jacobi matrix 


The study of the Sachs equation is simplified after performing a conformal transformation 
of the metric by a scale factor a, 

Spiv ® 9a u ■ (3.31) 

It can be checked that any null geodesic for g^ u , affinely parametrized by v , is also a null 
geodesic for g^, affinely parametrized by v with dt> = a 2 dv. The associated wave four- 
vectors then read A+ = a 2 k + Since the four-velocities of the comoving observers for both 
geometries are respectively u = d t and u = d^, so that S' 1 = a/, we deduce that 

co = g fll/ u tl k l/ = a~ 1 g fiu u fl k 1 ' = a _1 u). (3.32) 

The 3+1 decomposition of &+ is therefore 

k M = co (—vt 1 + n M ) (3.33) 

with n [l = a rffi implying = n^ja and 


ki 

fii = — . 
co 


(3.34) 


by 


The Sachs basis (n%) for the conformal geometry is then related to the original one (3.8) 


= ant 


h“=a 1 n". 


(3.35) 
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One can indeed check that the orthonormality (3.8) and the parallel transport conch 


tions (3.10) are preserved by the co nform al transformation with the use of the projection 


matrix S^ = a instead of Eq. (3.9). 


The separation four-vector ij^ is defined by comparing events only, independently from any 
metric. It is therefore invariant under conformal transformations. However, its projection 
over the Sachs basis changes (since the Sachs basis itself changes), indeed 


rf = = an a Jf — a fj a . 


n 1 n 

This implies that the Jacobi matrix transform as 

V nh = a V, 


ab • 


(3.36) 


(3.37) 


Hence, the angular distance Da in the universe described by a metric is just clDa, where 
Da is the angular distance in the universe described by the metric g /w . At lowest order in 
perturbations, k is the relative perturbation of angular diameter distance whatever is the 
metric used. As for the reduced shear, it remains unaffected by the conformal transfor¬ 
mation. In the remainder of this article, we will thus discard the effect of an overall scale 
factor, in order to simplify the computation. However it should be reminded that, as shown 
by Eq. (3.32), a conformal transformation has an effect on the energetic aspects of light 


propagation, that is on the relation between the redshift and the affine parameter U(v). 


3. General solution 


Using such a conformal transformation, it was shown in Ref. [68| that the Sachs equation 
can be solved analytically in a Bianchi I universe. This solution relies on the fact that the 
Sachs equation can be rewritten as 


d 2 V ab 

dv 2 




ac^cbi 


with the rescaled optical matrix given by 


(3.38) 


TZnh = W 


U! 


(O’ab)' + @ac&cb T ~ ®ab 


U 


The explicit solution of this equation is given in § VII.A-B of Ref. [68j . 


(3.39) 


IV. SMALL SHEAR LIMIT 


Definition 


The current observational status of the ACDM model shows that if the expansion is 
anisotropic, a/T~L has to be small. Moreover, since any primordial anisotropy is washed out 
by the expansion of the universe, the term C* in the evolution of the background shear is 


negligible compared to the integral term in Eq. (2.20). 


As discussed in the introduction, a late time anistropy may be generated during the 
acceleration of the universe, but the effect we are looking for needs to have an amplitude 
small enough to be below the detection threshold of ongoing observational surveys. 
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In full generality a linear and gauge-invariant perturbative expansion around an 
anisotropic background should be performed. It was developped in Ref. [22J [23] in the con¬ 
text of inflation and Appendix [A] derives the full perturbation theory for a post-inflationary 
era. While a numerical integration of these equations can be performed, it is clear from the 
previous arguments that an analytical insight in the regime a/T-L 1 is sufficient. 

We shall thus work in the small shear limit in which the background shear induced at 
late time by the anisotropic stress-energy tensor of the dark component is small, that is in 
the limit c/T/ <C 1. More precisely, we assume that 7 ^ — 5 ^ ~ 2 PiSij is a small dimen¬ 
sionless perturbation, and cr^/'H is of the same order as this homogeneous perturbation. 
We shall thus consider the Bianchi / spacetime as a homogeneous perturbation around 
an isotropic Friedmann-Lemaitre spacetime, hence ignoring non-linear corrections in the 
background shear as well. In order to implement this approximation scheme, we intro¬ 
duce a two-parameter perturbation scheme (see e.g. Ref. { 8 U]) in which, besides the usual 
Scalar-Vector-Tensor (SVT) perturbations over a flat Friedmann-Lemaitre background, the 
geometrical shear is considered as an extra perturbative degree. We refer to Ref. [ 8 T] for a 
detailed description of general Bianchi spaces in this approach. 


B. Spacetime description 

1. Metric 


We shall thus adopt the metric 


ds 2 = a 2 [—(1 + 2 $)d? 7 2 + 2B i dx t drj + ( 7 ij + hij) dx t dx j ] , 


(4.1) 


where h t j is defined as [see Eqs. (A19) and (A36) 

= -2 ( 7 ij + ^ + 2 Eij , 

and 7 ^ is here understood as the Euclidian metric plus a small perturbation 


crij da' oy, 


it K<1 


(4.2) 


(4.3) 


In order to simplify the notation, we also define the matrix 


Pij = diag(A) (4-4) 

such that 

7 ij = exp[2 \0\ij ~ Sij + 2, a { j = a) = ph . (4.5) 

Thus, controls the homogeneous perturbation. Indices are now raised and lowered with 
the Euclidian metric S t] and <W, and the vector modes Bi and tensor modes E t] satisfy 
d'Bi = d'Eij = E\ = 0. But since Pij is homogenous, everything happens as if we had usual 
cosmological perturbation, but also an infinite wavelength perturbations 2/3„ to the spatial 
metric. 

To control the perturbative series, we introduce the {n,p} notation, where n and p 
indicate powers in P and SVT variables, respectively. Thus, a term like Uij/l-i is of order 
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{ 1 , 0 }, terms like T and $ are of order { 0 , 1 }, while a product like a\d^/'H is of order 
{1,1}. However, since vector and tensor modes only appear due to the coupling between 
the shear and scalar modes [22], vector perturbations Bi and tensor perturbations E %] are 
also of order {1,1}. Hence, for any quantity A", one will consider the different quantities: 

• Arl°’°l ( 7 ^): the Friedmann-Lemaitre background value; 

• A^ 1,0 ^): the first order (homogeneous) scalar perturbed quantity in cr/7f; 

• A^ 0 , 1 ^(t 7 , x)\ the first order inhomogeneous perturbed quantity in T,...; 

• A't 1 , 1 }( 17 , x)\ the first order inhomogeneous perturbed quantity in both a/Ti and T,... 
and vector and tensor perturbations. 

Before moving on we should make some general remarks about the adopted perturbative 
scheme. Indeed, one might be worried that adding or (3^ as a small homogeneous 

perturbation to the background metric would not have any significant observable effect, since 
the SVT decomposition was already designed to describe the most general perturbation over 
a flat Friedmann-Lemaitre universe. Note however that SVT modes do not include a zero 
Fourier-mode in their spectrum (i.e., an infinite wavelength perturbation), since these modes 
will be isotropic by construction and hence merely rescale the background geometry. The 
tensor faj, on the other hand, is a homogeneous (i.e., space-independent) field, which by 
definition corresponds to an anisotropic zero mode. Thus, its effect cannot be absorbed in 
a simple rescaling of the scale factor. Moreover, this field sources the background dynamics 
through Einstein’s equations. 


2. Tetrad basis 


Given this expansion scheme, the tetrad basis associated to the perturbed metric up to 
order { 1 , 1 } is explicitely given by 


ey (*?-/y(i + *) + -d. 

0 s , ^ (6‘ + /%)(! - *) - d* + E ), 


H 


©o° = l-<k, 

0 a o = ! + $, 


0i° = O, 

& 0 = B\ 
Oo =~B\ 
©1 = 0 , 


(4.6) 


where d/ refers to the background spatial triad defined in Eq. (2.11). By chosing the 
observer to coincide with the timelike vector of the tetrad (tA = 0 O , u^ = — 0 -^) we obtain 
in general that the components of the direction vector n in the tetrad basis as 


k l - — -k%r , 


n l = QVG. 

P' 


(4.7) 


At the position of the observer, the direction of the geodesic in the tetrad basis n % 0 is also the 
direction in which the observation is made. Again, we remind that we are interested in the 
observables related to light propagation as expressed in function of this observed direction 
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Since we have introduced two types of tetrads ({t?} and {0}), there is an ambiguity 
whenever a tetrad index i appears on a tensorial quantity. First, for the geometric shear 
tensor, the tetrad index is defined with respect to the triad 'd p and we remind that a tJ = 


a lJ - = a 


j = diag(^') = /9h; see Eq. (2.12). Second, for partial derivatives the tetrad index 
corresponds also to the tetrad {$} and we define 


di = ~ d r - ftdj 


(4.8) 


ft makes clear the difference between a derivative in the direction of a tetrad vector di and 
the derivative in the direction of the vectors di associated with the Cartesian coordinates. 
Since the vector perturbations Bi and the tensor perturbations E t j are already of order 
{1,1}, there is absolutely no difference between their tetrad components Bi_ and Ejj at this 
order of perturbations and there is no need to be particularly careful. Everywhere else, a 


tetrad index refers to the tetrad {0} defined in Eqs. (4.6). 


C. Technical interlude 

Since we are interested in computing observables on the celestial sphere, spherical coordi¬ 
nates are much more convenient than Cartesian coordinates. This paragraph describes the 
use of such spherical coordinates in real space and of the associated derivatives (radial and 
on the unit sphere). Several definitions of covariant derivatives have to be distinguished. We 
finish by relating them to each other and to the spin-raising operator of spherical harmonics. 


1. Spherical coordinates in real space 


Consider a tensor depending on Cartesian coordinates Tq...*„(£*) (with indices raised and 
lowered respectively with and it can always be constructed by considering the tetrad 
components of a given tensor. In spherical coordinates, one can then define from the partial 
derivative di = d/dx 1 a covariant derivative D t on the unit sphere and a radial derivative 
d r . To be more precise, this requires the use of the projectors 


Sij — Sij XiXj , 



3 

with r 2 = y>*) 2 . 

i —1 


(4.9) 


Remind that S 13 = 5 tp S 3q S pq and Si = 5 ip S P j. The covariant derivative on the unit sphere 
•S' 2 (space of the Cartesian coordinates centered on the observer is denoted by Df 3 and is 
defined from the general projection 


- d r3 t • • 

r iy i 1 n-]n 


— gk gqi 


d 

CQn _ rji — TD 

■ °jn Q x k 1 1 - 9 " — r 


d_ 

dx> 


(4.10) 


where P[. ..] is to be understood as the projection of all free Cartesian indices with the 
projector Sj. This derivative only makes sense if the tensor itself is a projected tensor, that 
is, if it satisfies P[Tj 1 ...j n \ = T 

The radial derivative is then obtained simply by 


d 

g T ■ — x 1 —— 




(4.11) 
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Now, any combination of partial derivatives di applied to some tensor, can be decomposed 
in terms of radial derivatives d r and covariant derivatives on the sphere Df 3 . The simplest 
such decomposition is 


dif = x%f + -Dff , 


(4.12) 


for any scalar function /. The decompositions for projected tensors of various ranks is 


detailed in appendix D 1 For general tensors which are not necessarily projected, it is 


necessary to split them into their projected components on the sphere and their radial 
components before decomposing any derivative applying on them. Such decomposition for 
the vector and tensor modes is given in Eqs. (D5). To finish, it is easy to check that 


DfS jk = 0 , Dfe jk = 0 , 
where the completely antisymmetric tensor on the sphere is 

£ij = ^ijk% • 


(4.13) 

(4.14) 


2. Covariant derivative on the tangent space 


For any spatial tensor constant in space, such as cr y -, one can define scalar, vector and 
tensor holds on the unit sphere. First, one can define a scalar held on the unit sphere of 
observing directions, S' 2 | 0 b s , by contracting all free indices with the direction of observation, 


£ = -a i:k nln-. 


(4.15) 


Indeed, the observing direction can be considered as a point on S 2 | obs , whose spherical 


coordinates are (0 o ,<p o ), and £ from the expression (4.15) is thus a function of (9 0 ,ip 0 ), that 
is a scalar held on S' 2 | 0 bs- Then, to define a vector held on the unit sphere, one needs to 
contract one index with the observing direction and project the remaining one on the sphere. 
Furthermore, in order to get a tensor held on the unit sphere, we shall project the two free 
indices on the sphere. These projections are obtained by contraction with the screen basis 
vectors n° at the observer. For instance, the vector and tensor helds on the sphere build 
from the geometric shear are simply 


£ a = n° J a r jil £ a & = n^-nl'-a^ . 


(4.16) 


We remark that £ a b, which is a symmetric 2x2 matrix, is not traceless. In fact, using 
the partition of the identity W = n-n/ + n\ n\ + n\n^ the trace is given by = 

—cTikn-n- = —2£. Alternatively, the vector and tensor helds (4.16) can be obtained by 
applying successively the covariant derivative on the unit sphere D a to £( 6 ) 0 ,(p 0 ). Indeed, 
with this method, we find the relations 


S« — DqYj , 

Sab = D a DbTj + 2 habS , 

(4.17) 

D a D a T, = — 6 £ , 

S(ab) ■ 

Note that the metric and the antisymmetric 

tensor on the sphere are obtained from 


Sab = 2 

e a b = 2m|j - np 

(4.18) 

and satisfy 

D a $bc 0 , 

D a ^bc 0 . 

(4.19) 
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3. Background geodesics and identification of covariant derivatives 


The covariant derivative D a , related to the unit sphere in the observer’s tangent space, and 
the derivative Df 3 , related to the unit sphere of Cartesian coordinates, are fundamentally 
different. But, they can be related in a simple way. Indeed, the solution to the background 
geodesic at order {0,0}, that is, the geodesics of the spatially flat Friedmann-Lemaitre 
spacetime, is given by 


n- = n°- 

IL a n a 


dx*{0’°} 


x i{m = X <, 


(4.20) 


where, we remind the reader, nf is the direction of the geodesic at the position of the observer 
in the tetrad basis. This is the direction of observation, since we have oriented the geodesic 
toward the past. 

There is thus a straightforward identification between the sphere of the directions of 
observation, lying in the tangent space at the observer (the set of directions S 2 \ 0 b s spanned 
by n^), and the set of points of R 3 reached at an affine parameter x (or v) on the background 
geodesic. Indeed, the points spanned by the coordinates a t a given affine parameter 

X are such that 

^(CWViO.O) = r 2 (x) = x 2 

and form a sphere in the Cartesian coordinates. We can then subsequently identify this 
sphere of radius x to the unit sphere S' 2 | space . 

This means that we can identify n\ with x 1 and then D a on S' 2 | 0 b s with 7i°-Df 3 , the 
projection onto the screen basis n° being used only to switch from the extrinsic point of 
view of the derivative (the projection of the Cartesian derivative onto the sphere) to an 
intrinsic point of view on the sphere. In the rest of this article we thus replace the notation 
Df by A, nfDf by D a , and n° ± L Df by D±. 


4■ Link with spin-raising operator and spin-weighted spherical harmonics 


The covariant derivative on the unit sphere is related to the usual spin-raising and spin¬ 
lowering operators. In spherical coordinates, these operators are defined for a spin-s quantity 
by 


0X S 

0X S 


- siiT 9 

- sin" 5 9 


1 

sind 
1 

sin 9 


d e + i—— 


de - i—— 


(sin ~ S 9X S ) , 
siir s 9X s ) . 


(4.21a) 

(4.21b) 


They are related to the covariant derivative through 


= -V2n°_ a D a = -y/2 D+, 


= -V2n° + a D a = -y/2 , 


(4.22) 


the vector n° ± being defined in Eq. (3.11). Hence, for a tensor field of spin +\s\ on the 

Vl...V a _ 

-V2X^X l 


sphere, X Vl - Va = X s n° + Ul ... n° + v % and a tensor of spin —|s|, Z Vl " Ms = Z s n°_ Ul ... n°_ v % we 
have 


nYn° + Ul .. 

• <H( 

0X s )n°Jn° + Ul ... n° + Vs 

(4.23a) 

, )n°/n°_ i ' 1 . 

.. n°_ Vs + 

(0Z- s )n°Jn OVl ...n°_ Va ■ 

(4.23b) 
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Since the spin-weighted spherical harmonics satisfy the property 


YL = 


(^ + s)! 


fY p 


£m 


if 0 < s < 


-l) f 


(l + s)! 


(4.24) 


Y, 


£m 


if 


< s < 0, 


(t-s)'- 

any number of covariant derivatives applied on a spherical harmonic can be computed using 
the properties (4.23). 

As an application, consider the expansion of the variable E in spherical harmonics 

+2 

E(xX) = E m(x)Y 2 m(ni ). (4.25) 

m =—2 

If we align the azimuthal direction with an eigendirection of the geometric shear, the mul¬ 
tipoles coefficients are then given by 

S 2o(x) = ~\J\\PM + &M\, S 2±2 (x) = \^\PM - PM\ • (4-26) 

The most useful derivatives are then easily obtained to be 

£±(xX) = n°-n Q ^(x) = AfcS(xX) = tV% E ^ 2 m(x)Y^(K ), (4.27a) 

m 

E±±(xX) = = D ± D ± Y(x,<) = V 6 ^ 2 M)Y 2 ^(<) ■ (4.27b) 


Similarly, if we expand a scalar held </?(%, x®) in spherical harmonics 

vixi x% ) = E^(x,r)r^(i*), 

£,m 

then the most useful derivatives are 


At¥>(x,z*) = E^m(x,r)y^ 

£,m 


±1 






AfcAfc¥>(X> ®*) = o 


1 /(£ + 2)! 


2 I/ (£ — 2)! 






£,m 


(4.28) 


(4.29a) 


(4.29b) 


2 j D ±j D tV ?(x, = D a D a <p( X , j) = + 1) E ^(X, r)y, m (^). (4.29c) 

£,m 


D. Geodesics and Sachs equations in term of the Friedmannian coordinates 

In the approximation that we are considering, we can solve the perturbation equations 
and the Sachs equation up to order {1,1}. We shall define the distance down to the lightcone 
on the Friedmann-Lemaitre background spacetime as 

X = Vo ~ V■ ( 4 -30) 
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The geodesic equation (3.1) takes the form 


d k v 
dv 


+ T'' aB k a k f} = 0, A;° = -^, 


du 


and using (4.7) can be rewritten directly in terms of tetrad components as 

d k- d ki 

dv dn 

d k°- _ _c ikg 
dv 


“TT = “TT = (k-) 2 (uJkjinhtl + - u M n l - uJ m n J - 


= (k-) 2 [ -iv m nhil + u m n- 


du 


(4.31) 

(4.32a) 

(4.32b) 


where the affine connections are defined in Appendix B 1 Instead of the parameter v we 
shall use the parameter x since, once the wave-vector is integrated, we have 

% = (433) 
The position on the geodesic then becomes a function of the parameter x and the initial 
direction n- 0 . Finally, the Sachs equation with the parameter x reads 

1 


d 2 V ab 1 dk°dV ab = _ 

dy 2 k° dx dx ( k°) 


arV , 


ac^cb • 


(4.34) 


V. ANGULAR MULTIPOLE CORRELATIONS IN ANISOTROPIC SPACES 


As previously explained, we adopt an observer point of view in which all observable 
quantities are considered as functions of the direction of observation n Q and of the affine 
parameter v or, equivalently, of the redshift z, keeping in mind that the later also depends on 
n Q . All these quantities can be decomposed on a basis of spin-weighted spherical harmonics, 
Yf m . The goal of this section is to derive a set of formal expressions concerning these 
expansions and to establish general results of the two-point correlation function valid in 
Bianchi / geometries. 

We consider that the universe has undergone an early period of isotropic expansion fol¬ 
lowed by a late-time anisotropic phase. This is in sharp contrast with the approach of 
Ref. [30] . in which the universe is supposed to have an early inflationary stage followed by 
an isotropic evolution (so that geodesics are Friedmann-Lemaitre geodesics and anisotropy 
is imprinted only in the source term). 

The tools we shall develop are not specific to weak-lensing and can be used in other 
contextz, such as the study of the cosmic microwave background. We first describe, in 
§ V A the general expansion of spin s quantities. This will allow us to express their angular 


power spectrum in § | V B We conclude by demonstrating that, while spatial parity symmetry 
implies that the EB correlation matrix vanishes, some off-diagonal correlaions are necessarily 
non-vanishing and encode information on the geometrical shear. 


A. Multipolar expansions 

The spin-s polarization basis is defined as a tensor product of s polarization vectors as 


s _ . 


m G = 


n, 


n2 


n, 


if s > 0 , 
if s < 0 . 


(5.1) 















Under the action of an active rotation R, this basis transforms as 

R[m s 0 \ = R ■ m^R- 1 ■ n Q ), (5.2) 

where n Q is the vector along the line of sight at the point of observation. Spin-weighted 
spherical harmonics transform under the same rotation as 

R[Y£Jn 0 )m.i\ = YfJR- 1 • n,)«K] = V Y^,(„ 0 )rn‘D‘ m ,JR ), (5-3) 

ml 

where are the components of the Wigner D-matrix. This means that they trans¬ 

form like normal spherical harmonics provided they are accompanied by the polarization 
basis to which they are associated. 


Now, any cosmological observable X s of spin s can be expressed in the form 

'■Xs 


X'(x s ,n„ )m‘ 0 = I S x ‘(x a ,X,n.„)ml dx, 


(5.4) 


where x s refers to the position of the source. Note that we are explicitly making use of the 
small shear expansion, since the source term is integrated along a geodesic of the Friedmann- 
Lemaitre spacetime. This means that in order to compute A" s at order {n,p} one needs to 
determine the source S Xs at the same order. The source term S Xs (x s ,X, n o) has to be 
understood as 

S Xs (x a ,X,n 0 ) = S Xs (x s ,X,x\n 0 )\ xi=xn i , (5.5) 


that is, evaluated on the background geodesic. Moreover, thanks to Eq. (4.30), the parame¬ 
ters x and y s can be thought as time coordinates. The intrinsic angular dependence of S Xs 
on n 0 is a consequence of the (possible) non scalar nature of the source. Moving forward, it 
is convenient to decompose X s into spherical harmonics as 


X s ( Xs ,n 0 )m s 0 = J2 x L(x s )Yln(no)m s 0 , 

£,m 


(5.6) 


which will allows us to define multipolar correlations at unequal times of the form 
(A"! m (x S i)A^'m'C\s 2 ))- I 11 or der to compute these angular correlators, we first need to take 
the Fourier transform of the source (5.4) off the line of sight 


S' s {x s ,X,x\n 0 )m s 0 = 


d 3 k 


S a {x s ,X,k,n 0 )e 


ikx 




(5.7) 


(2jt) 3 / 2 

in the sense that we do not bind x l to x by the relation x l = x n o an d X das to be thought 
as a time coordinate thanks to Eq. (4.30). Then, the intrinsic dependence of the source on 
n D is further expanded in terms of spherical harmonics, with the latter being defined with 
respect to an axis aligned with the Fourier mode k. That is 


S Xs (x s , X, , k, n 0 )m s 0 = ^ S£(x s , X, k)i e ] j^^R k [Y; m (n 0 )m s 0 \ , (5.8) 

i,m 

where R k is a rotation that transports the azimuthal direction to the direction of the Fourier 
mode k (see Appendix [C| for details about this notation). The terms with m = 0,1,2 
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correspond here to scalar, vector and tensor perturbations, respectively. If we now make use 
of the Rayleigh expansion 


e ikx = 4vr Sjt(kr)Y; m (k)Yz m (n 0 ) = i e y/(4n)(2£ + 1 )j e (kr)R k [Y m (n 0 )], (5.9) 

t,m l 


with r = SijX l x J , and insert the decomposition (5.8) into Eq. (5.7), we find, after comparing 
Eqs. (5.6) and (5.7), that 


XL(Xs) = \- 


d 3 k r AX '^D! M > 

^ m I 


■e A fif. + 1) Y^ s .(«) /U nx, 


47T 


Y ‘it ” \ kr ) s im’{x S ,X, k) 

(5.10) 


where we have introduced the definitions 


s .{t'm')( X s^m'Oin' ■ / \-L+l'-t ( 47 r)( 2 L + 1 ) 

3 , W-y jiWi Y(2C+ 1 )(2f + l)’ 


(5.11) 


and 


s «a m3 = J d2 OC 1 (n 0 )R 2m2 (n 0 )y/ 3ffl3 (n 0 ). (5.12) 

The dynamical evolution and the initial conditions of the source can be split as 


St(Xs:X,k) = T^(x s ,x,k)$i(k) , (5.13) 

where is the (anisotropic) transfer function and <hj(fc) is the primordial gravitational 
potential. Then, assuming that anisotropies are induced at late-time evolution only, the 
statistics of the primordial power spectrum must obey 


(^(fc^ffa)) = P(k)S 3 (k - q ), 


(5.14) 


with P(k) being the (isotropic) primordial power spectrum. To account for the angular 
dependence of the transfer functions, we further decompose them as 

Tt(x„X,,k) = Y ^“'(Xs, X, k)Y LM (k ). (5.15) 

L,M 


B. Expression of the two-point angular correlators 


These formulas can now be combined (using in particular Eq. (D25) to integrate out all 
spherical harmonics) to give an expression for the correlation between the multipoles of two 
different observables A" Sl and Z S2 . We find 

2 r°° 

W,' mi (x s .)ZS;,(Xs a )) = -(i)' , (-i) 4 / cl kk 2 P(k) 


7r 


dxi dX 2 £E 5] (-ir^"'AY ,} (kxi) 

e,m V x ,m' x L u Mi 
l' 2 ,m' 2 l 2,m 2 

X -T, L t 'rX 1 (Xs,’Xuk) z -^7*(x». X2, k ). (5.16a) 
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A central quantity in this description is the two-point correlation function of the E- and 
5-modes of a given spin-2 observable (as, for instance, the cosmic shear 7 ). This expression 
requires the decomposition of A" ±s , ±s jf ^ and Sf^ in their even/odd parity pieces as 


x L(x) = E^{x) + isgn(s)B-^{ X ), 

±s jf m \x) = |s| ef m \x) + isgn(s) |s| / df m \x ), 


s hn(x a ,X,k) = l T fm 3 (x S iX,k) + isgn(s)T^;' s (x s ,x,fc)]$ i (fc). 


(5.17a) 

(5.17b) 

(5.17c) 


Note that a spin s = 2 field will have both E- and 5-modes, while a scalar (s = 0) field will 
only have the E mode, so that Q fif m ] = 0 and Bf n ° = 0. 


From these expressions and Eq. (|5.10|), one can verify that 

Efa(Xs) = 


d 3 k f 

J 0 / 

m ' 


(2 i + 1 ) 


47T 


(5.18a) 


B lm(X s ) = 


E [‘'€ m,) (kx)T^(x „,*,*) - m ‘\kx)T^(x„x,k) $,(*), 

/d 3 fc f 3 dxEDL,ARk)i ,J{2e + 1) 


(5.18b) 


E (Xs,X,k) + l 7<f’”' ) (WTS'fc.X.fc)l Ki(fe) ■ 


Then, the EE and 55 covariance matrices can be computed by simply taking appropriate 
combinations of XfA- In order to simplify the notation we define 


M 


ACA'C' 


(Xsi 1 XS 2 ) 


9 „ r°° rx si rx S2 

— (i) £l (—i)^ 2 / d kk 2 P(k) / dxi / dy 2 (5.19a) 

E E E (-1 r i+mi A '^ mi> (k Xl ) 

t,m Li,Mi 

£' 2 ,m ! 2 l 2 ,M 2 


-m\ rxm\mM\ 


”AxAT.;s<Xsi>* i>*) 


C1J-1L2M2,* 

1 e’ 2 m ' 2 


(Xs2 j X 21 k') 


together with the notation 


AAl'rri) 

■h 


Wef n,) if A = E Xs , 

H/f ' m>) if A = B x *. 


Thus, the EE and 55 correlations become 


(E, 


X s rpX s * 
t\mi(. 2 m 2 


> = E M f AC 


m\l 2 m 2 


A,C 


(B. 


X s rx V s* 
<imi ^£ 21 x 12 


\ = v ka acAC 

1 / j £\m\l 2 m 2 


A,C 


IfiSE + , 

[&A&C + + <^] , 


(5.20a) 

(5.20b) 


where, in the last equality, we have introduced the notation according to which A equals 
E Xs (resp. B Xs ) whenever A is equal to B Xs (resp. E Xa )] the same holding for C. The EB 
correlation can be computed using the same method. 
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The expression (5.20) is quite general. Let us first focus on its diagonal part, which can 
be characterized by the following estimators of 


Cf E (x„,xJ = 2 ^ Estate,)£&*(*»)>, (5.21a) 

m 

cr(x 3 „xj = ^£<B£(x sl )B£*(x 91 )}. (5.21b) 


The angular power spectra are then given by simpler expressions 


C e (X S 1 ,X S2 ) = - 


x 


x 


d k k 2 

47T 

E 

. h ' 

E 

. ii ' 


p(k) Yl \- 6 * s c+tit 

m,L,M,A,C 


B xE xB 
C — b A bc 


^A^c] 


"Xsi 


d Xi A TiZ(x si, X S2 , k) A jt xm> {kx i) 


A-(tim), 


^*Xs 2 


dX2 C TS( Xsi , Xs2 ,A;)^r m; fe) 


C-(fem) , 


(5.22a) 


with, again, a similar expression for the 5-modes. 

The case where the transfer functions are isotropic is easily recovered. To see that, let us 
consider the simpler situation where A " S1 = Z s2 = 0, with 0 being the CMB temperature 
fluctuations. Since 0 is a scalar, then it is a pure 5-mode with no 5-mode. In the previous 
formalism, we just need to set E T and B T[Z = 0. Then, using e T^ f = 
V^T° m 6 L0 5 MO , we get from Eq. (|5.16[) the standard result 


Cf 0 = - 

7r 


k 2 dkp(k ) ^ 


d x' T ®rn( k , y!)jf m) C kx' 


(5.23) 


where in this specific case it must be understood that the visibility function is included in 
the transfer functions TZ(k,x')- If we consider only scalar sources, then only the m — 0 
mode contributes. Analogously, if we also have sources with no intrinsic direction (like, for 
example, no Doppler effect in the CMB), then we have i' = 0. 


C. Implication of spatial parity 

We would like to briefly elucidate the relationship between the symmetries of the under¬ 
lying background spacetime and the cross-correlation functions of different observables. In 
particular, we want to show that (spatial) parity symmetry implies that the diagonal piece 
of the EB correlation matrix is zero, while off-diagonal terms may not necessarily be. 

We start by noticing that under a parity inversion {x, y, z} -A- {—x, —y, —z}, or, equiva¬ 
lently, {n o ,n 0 ,n^} —$■ {—n Q , — the polarization vectors transform as 

n o( n o) -t ~ n o ( n o) ■ (5.24) 

This implies that the polarization basis should transform under parity as 

m s 0 {n 0 ) -x- (-1 ) s m~ s {n 0 ). (5.25) 
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Moreover, the sources transform as 


S a {x S iX,k,n 0 )m s 0 (n 0 ) ->■ (-1)*S s (x s ,X, ~k, -n 0 )m 0 s (n 0 ). 


(5.26) 


We now demand that any physical quantity remains invariant under a full parity inversion. 
That is, if at the same time we transform k —> —k, n Q —» — n Q and m s 0 —> (—1 ) s m~ s , 
the source S Xa of a physical observable X s should remain invariant, which from Eqs. (5.7) 
and (5.26) implies the condition 


s •(x a ,X,k,n 0 ) = (-1 )’S -(x s ,x,-k,-n 0 ). 


(5.27) 


If we take into account the parity transformations of the Wigner matrices and of the spherical 
harmonics - see Eqs.(D9 Dll) - then a comparison of the previous expression with Eq. (5.8) 
shows that the parity condition translates to 

St(Xs,X,k) = (-1 r +s S?~^ Xs ,X,-k). 

Then, we rewrite Eq. (5.10[) as 


(5.28) 


X’,jx s) = \jlj d 3 kx; m ( Xs ,k)Hk). 


(5.29) 


The above expression should be seen as a definition of Xj! m (x s , k), and corresponds to the 
contribution of each Fourier mode to the observable, but it is not its Fourier component. Its 
expression can be obtained by plugging Eq. (5.13) into Eq. ( |5.10[) a nd then comparing with 
Eq. (5.29). Then, if we impose the symmetry (5.28) to Eq. (5.29), using again the parity 
transformation of the Wigner matrices, we finally find that 

X’ tm (x s, k) = (-l)‘ + ‘X^( Xa , -k) . (5.30) 


From this expression, it is straightforward to check that the E- and 5-modes of a spin-2 
quantity transform under parity as 


e£(x, k) = 1 (X]Jx, k) + x£( x , k)) = (-1 )‘e£( x ,- k) (5.31a) 

B£(x,k) = Y(XUx,k)-X£( X ,k)) = (-l) ,+ ‘B£(x,-k). (5.31b) 


We can now establish our main conclusion: given the above symmetry, together with trans¬ 
lational invariance of primordial fluctuations (see Eq. (5.14)), it follows that, for the E and 
B modes of a spin-2 held, we have 


Kl»(x)B?X,(x)) = l J d 3 k (E£ m (x,k)B%; 2 (x,k))p(k) 

= (-1 )' ,+fe+1 l J d 3 k (E-£ m (x, -k)B%Xx, ~k))P(k ) 


= (-1 ) fl+fe+ D J d 3 k(E£ mi ( x ,k)B£ m \(x,k))P(k) 

= (-l ) f ‘ +fe+1 (Efc, t (x)Bg£(x )), 


(5.32) 
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where, from the second to the third line, we have used d 3 fc = / +c ^d 3 (— k). Similarly, 
one can show that the EE and BB covariance matrices obey 

= (-1 ) h+h (E£ m E*Q, (5.33a) 

= (-1 (5.33b) 

We have thus proved that correlations between E and B modes will vanish whenever £ 1+^2 
is an even number. In particular, the diagonal part of the EB covariance matrix is always 
zero in spacetimes that respect parity (but not necesserily isotropy). Evidently, the same 
holds for the multipolar coefficients of a spin-0 quantity, such as the CMB temperature a ^ m s 
(see e.g. Ref. [82]). 


VI. PERTURBATION SCHEME IN THE SMALL SHEAR LIMIT 


A. Expansion scheme 


The structure of the computation has been detailed in 
order, we need to 



Let us recall that order by 


1 . solve the geodesic equation perturbatively in order to determine the displacement from 

the reference Friedmann-Lemaitre geodesic, n°), and the local direction of 

propagation, n^ n ’^(x, n°). Note that x l ^ n,p ^(x, n°) is split in a radial component, 
hr l ' n,p ' ! (Xi n°), and an orthoradial component which will be related to the deflection 
angle a a ^- n,p i (y, n°); 

2. determine the transport of the Sachs basis, ni n ’^(x,n°); 

3. expand the Sachs equation and determine the source terms for 'D f |j)’ J) ’’(y, n°); 

4. determine the evolution of the perturbations at the required order; 

5. perform the multipolar expansion in terms of the direction of observation n°. 

To avoid confusion, we shall use the notation that X^ n,p ^ includes all terms up to order 
{n,p} while 5X^ n,p ^ contains only the terms of order {n,p}. 


B. Order {0,0} 


Since on the background (i.e., a Friedmann-Lemaitre spacetime), the metric is just the 
Minkowski metric, thanks to the overall conformal transformation described in § |III D 2 


the (conformal) Riemann tensor vanishes, so that 7^’°^ = 0. Since the wave-vector is 
decomposed in accordance to Eq. (3.2), in which we can always choose to set U D = 1, one 
deduces that it is given by 

jfeO{o.o} = _ 1 , k i{Qfi} = n i ( 6 T) 


The Sachs equation (4.34) trivially reduces to 


d 2 £> { ?’ 0} 

-^— = 0 

dx 2 


( 6 . 2 ) 
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so that the Jacobi matrix is given by 


PS 101 = Df°\ X )U, D[ m ( X )=X, (6.3) 

and its components reduce to 

k {O,O} =7 {O,°} = v{ O,O} = 0 ' ( 6 . 4 ) 

This completely specihes the property of the geodesic bundle at the background level. 


C. Order {1,0} 


At this order, the spacetime remains homogeneous, but it now has an anisotropic pertur¬ 
bation described by the shear a^, from which we can define a scalar held E on the 2-sphere 
by 

S C\) = ^ik(x)K n o- 

We also introduce a new scalar function, 

E(x) = ^iMK n o, 


where f3ij is dehned in Eqs. (4.4) and (4.5). 




(6.5) 



(6.6) 

in terms of 

these two 

in § 

IV C 2 

In what 


follows we shall use the convention B 0 = B(x = 0). 


1. Geodesic equation: tangent vector 

At this order, the 4-velocity of a fundamental observer is just u^ = (d^) M = , so that 

(6.7) 


U = k % - k° - 


From Eq. (|4.3l|), and using the fact that, at first order, the only non-vanishing Christoffel 

<srg 1>0} = a) , (6.8) 

dk g{ifi} dk o{i,o} 


symbols are 
we obtain that 


xn 0{l,0} 

= on > 


dy 


dx 


= -Oiknln* 


It thus follows that 


fcSli’O} = A: 0{1 ’ 0} = -1 + 2 [B(x) - B 0 \ = -1 + Jfc 0{1 


,0} 


(6.9) 


( 6 . 10 ) 



events of emission and reception. 
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The evolution of the spatial components of the wave-vector is easily obtained using the 
tetrad components first. From Eq. (4.32) we obtain 


k- {l,0} (x) = < + [ v-Ax'hodx' = n- 0 + [^.(0) - P)(x)\ n o (6.11a) 
Jo 

k i{1 ’ 0} (x) = ^ {1 ’ 0} (y) - P)(x)ni = n% + [P)( 0) - 2 /^(y)] r& (6.11b) 
kl m (x) = + Pij(x)n o = < + P l j(0)rio (6.11c) 

Again, this corresponds to the small shear limit of our previous general result (68] . where 
ki = ki(0 ) is used first. 


2. Geodesic equation: real space 


The parametric equation of the geodesic is obtained from the integration of Eq. (4.33) 
at order {1, 0}. Separating the difference between the position at order {1, 0} and the posi¬ 
tion of the background geodesic into a radial displacement and an orthoradial displacement 
according to 


x^ 1,0 ^(x) = x n l + foe®* 1,0 * with = nl5r^ 1,0 ^ + ya^ 1 ’ 0 -^®, (6.12) 

which defines the deflection angle aA^ 1,0 }, and where we have used that x^°’°^(y) = yn^, we 
get 

^ {1,0} (y) = -2 [ B(x r )d x , (6.13a) 

Jo 

« a{1 ’ 0} (x) = D a B 0 - - P (6.13b) 

X Jo 

in which D a B 0 stands for (D a B ) x=0 , and where the last equalities of the equations above 
made use of an integration by parts. Note that 


lim-= n- - Pij(0)ni ■ (6.14) 

x->o y JV 


This can be interpreted simply, because it means that very close to the observer, everything 
happens as if d l3 is constant and equal to /%(0). Thus a constant change of coordinates 
x 1 = x l + /3*(0)x J transforms the metric from X ij = &ij + ‘jPijid) to the Euclidian metric Sij. 
The geodesic in these new coordinates is simply the Euclidian one, x*(y) = X n o- 


3. Evolution of the direction and screen vectors 

The infinitesimal change of a unit vector lies in the plan orthogonal to it. The perturba¬ 
tion of the direction vector is thus of the form 

7T {1 ’0} = < + ® a{1 ' 01 n"'. (6.15) 
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From the previous results for k- and k-, we get immediately 


n 


= n\ + w a{m n°j -, cu a{1 ’ 0} = / D a E dy'. 


(6.16) 


The transport equation for the screen basis is in turn given by 


i{l,0} _ oi _ i { 1,0} 
IL a IL a IL o ^ a i 


(6.17) 


and it can be checked that the screen basis (3.10) does remain orthogonal to the direction 
vector. 


4- Sachs equation 

At order {1,0}, given that 77^' 0/ = 0, the right hand side of Eq. (4.34) reduces to 

= K M xSte = xK im , 


so that the Sachs equation (4.34) reduces to 


d%>V° l _ d»“C'°> {m 

-d ab + xl<- ab 


dy 2 


Its first integral yields 


d5V 


{1,0} 

ab 


dy 


p * fddk 0 ^ 


(6.18) 


d* =s ^k 

The first term gives 5 a b[l + 5fc°^ 1 ’°J'(y)] = S ab [2 + fc°^ 1,0 ^(y)], so that 

sv ll ,0} (x,n°) = Sab [ [2 + fc 0{1 ’ 0} (y)] dy + [ dy '[ y^’° } (y , n°)dy . (6.20) 

Jo Jo Jo 

The double integral on the right hand side can be performed by means of an integration by 
parts. This gives 

x 


(6.19) 


dy' / y^Ji’ 0} (y,n°)dy= / y(y - y)^’ 0} (y, n°)dy 


'0 


rx rx 

/ dX l 

Jo Jo 

from where we finally conclude that 

rx 
10 

As detailed in Appendix |B 2 \ the source term takes the form 

1 

2 (°*j 


sv ii 0} (x, n°) = r { [2 + k°W(x)} S ab + (y - y) y^’ 0} (y, n°)} dy . ( 6 . 21 ) 


^S’ 0} (x, n°) = -^«)^ + ngn^y)' = -<UX' + D {a D h) V . 


( 6 . 22 ) 
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By inserting the above in Eq. (6.21) and using Eq. (6.11), we find the following expressions 
for the convergence and shear, 


k {1 ’ 0 } (x, n°) = - f X [2E + *£'] dy , (6.23a) 

Jo X 

7 £’° } (x,n°) = f X {X ~ X)X D {a D b) E'dx, (6.23b) 

Jo X 

which simplify to 

K {1 ’ 0} (y,n°) = -B(x) ~SB 0 + - [ B(x)d%, (6.24a) 

X Jo 

72’° } C\, n°) = D {a D b) B 0 + D {a D b) B - - P D {a D b) B{xW • (6.24b) 

X Jo 


Note that in the limit y —>■ 0, —> 0 and 7^’°^ —> 0 as it should be, given the initial 

condition (3.16) for the Jacobi matrix. Finally, since the Jacobi matrix (6.20) is symmetric, 


at this order we have 


U {1 ’ 0} = 0. 


(6.25) 


We have checked that these results match those found in (68] when expanded in the small 
shear limit (where the special choice B 0 = 0 is made). 

We can now perform the expansion of these observable quantities in terms of spin- 
weighted spherical harmonics. Using the results of § |IV C 4( the shear can be projected 


into the hclicity basis so as to transform the covariant derivatives into spin-raising and spin¬ 
lowering operators. The spherical harmonics components of the convergence k and of the 
cosmic shear 7 ± are then easily obtained at order {1, 0}. Both reduce to a quadrupolar con¬ 
tribution, inherited from the quadrupolar contribution of E, so that their only non-vanishing 
coefficients are 


, {1 > 0} 

^2m 


ix) = -3fl 


O 

2 m 


£>2m H- 


X Jo 


#2m(x)dx, 


7?i 1,0} (x) = V6 ^2m + #2 B 2m{x)&X 

We conclude that = 0 and 


(6.26a) 

(6.26b) 


D. Order {0,1} 

We follow the same method for the order {0,1} as for the order {1,0}. This corresponds 
to the standard approach to weak-lensing in the linear regime of cosmological perturbations. 
Our main goal is to rederive these standard results in our formalism, so as to serve as a 
basis for the study at order {1,1}. Note that, at this order, we only need to include scalar 
perturbations since, as stressed before, vectors and tensors modes are of order {1,1}. 
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1 . Geodesic equation: tangent vector 


Using the definition of the deflecting potential as p = $ + T, the energy of a photon 
evolves according to 


dfc-i 0,1 } dip d<f> d£0{o,i} dp 

dy dx dx ’ dx dx dy ’ 

the solution of which is 

£0(0,1} = £0(0,1} + $ = _i + 2$ - - [ X ^dy = -1 + 5k°W , 

Jo dx 


(6.27) 


(6.28) 


where it is understood the integrand is evaluated on the background geodesic, i.e. that 
x l = x n o> a thne associated with y ■ The notation is intentionally simplified in this 
section, so, for example, fc°d°’ 1 l means /c°^ 0,1 ^(x, n^), $ means <h(x,a;*) with x l = x n o-> an d 
so on. In other words, it is understood that everything is evaluated on the background 
geodesic at parameter y. The only exception is <h 0 , which is the potential <h evaluated at 
the observer, that is, at y = 0. Note that the total derivative d/dy, i.e. the total derivative 
along the background geodesic, satisfies 


since, from Eq. (4.30), 


dp dp . _ „ 

— = — + p = -p + x dip , 

dy dx 


dp__, 

dx ^ 


(6.29) 


(6.30) 


2. Evolution of the direction and screen vectors 


The spatial component of the vector k M evolves according to 


dfcd°’i} ^ . dT 

—-= -dip + n-—, 

dy dy 


dk i{ °d} dT 

—-= -dip + 2n— , 

dy dy 


the solution of which is 

£*{o,i} = £*{o,i} + =n i [i + 2T(y) - T 0 ] 


dipdx • 


Using Eq. (4.7), we then deduce the evolution of the direction vector 

dn-(°’i} 


dy 


the solution of which is 


= -S^-djp, 


’ } (y, <) = < + ^ 0>1 >(x, <)<> t j7°i°> 1 >(y, n i) = 

Similarly, the evolution of the screen projectors leads to 


D a pdx • 


n- 


*{ 0 , 1 }/ 


(x, K) = C - ™ a{ 0 , 1 } (x, 


(6.31) 

(6.32) 

(6.33a) 

(6.34a) 

(6.35) 
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3. Geodesic equation: real space 


We can then determine a;do,i} from Eq. (4.33) using 

'dx®\ ^°’ 1 ^ 


d XJ 


k 


= *(A) = 4(1 + v)-J SVdjipdx 


{ 0 , 1 } 


and this leads to 


x 


do,i} = xn i + n l - [ (pdx ~ n° a ~ [ d ^ X . ^ D a (p . 


' 0 


'0 


X 


(6.36) 


(6.37) 


4■ Sachs equation 


Finally, at order {0,1}, the right hand side of Eq. (4.34) is simply = x7?Ic 0 ^. 

Thus, the Sachs equations becomes 

^2 - dx ° ab + ^ ,l ab • ( b -d») 

The solution of the Sachs equation follows formally the same steps as in the case {1,0}. 
That is, it can be integrated twice, and after an integration by parts for the double integral 
over the Riemann term, we get 

M> 2 ’ 1 } (x,n°) = I"" { [2 + fc'X°’ I >(S)] S ab +(x~x) xKJHx, n°)} d* . ( 6 . 39 ) 


Now, using the perturbed expression for JZ ab found in the Appendix B2 (with = 0), 

(6.40a) 


X 2 K% 1} = x 2 n°a n b L [-didjip - % (T" - 2 n-d^' + rEn^d p d q ^)] 


= -D {a D b) ip - 5 ab x 


2 ^did l (p + T" — 2n-<9;\E + ^ n p n q d p d q {T — <f>) 


(6.40b) 


and the expression d/c 0 ^ 1,0 ^/dy given in Eq. (6.28), one obtains the formal solution of the 
Sachs equation (4.34) as |62 

i rx i rx 


5V 


{ 0 , 1 } 

ab 


= x 


8ab[l-V(x)-®o+- <p(x) dx-x 


X Jo 


x-x 

lo XX 


D c D c v{x) dx 


X-X 


D( a D b) ip(x)dx 


(6.41) 


Jo XX 

from which k and 7 a b can be read directly from the expression in brackets in the first and 
second lines respectively; see our definitions in Eq. (3.17). Note that, since there is no 
antisymmetric part in we conclude that l/f 0,1 ! = 0 . 

Dropping the (unobservable) monopole correction due to the local potential 4> 0 , we get 
their multipoles as 


4m 1} = -dbm(x) + - / 7km(x)dx + 

X Jo 


^ +1) f x X~X 


lo XX 


±{°,!} 

i£m 


1 {£ + 2)! f x X~X 

' 2 V (^-2 y.Jo XX 


<Pim(x) dx- 


<P£m(x) dx, (6.42a) 


(6.42b) 


From which we conclude that = 0 and = y^ 1,0 ^ 


tm 


Her, 


40 






















5. Angular power spectra 


To determine the angular power spectrum of the convergence k and of the .E-modes of 
the cosmic shear, we follow the procedure described in § V At the order {0,1}, the transfer 
function is isotropic and there are only scalar sources. Consequently, only E-modes are 
generated. The power spectrum for the E-modes is then just given by 


c, 


EE{ 0 , 1 } 


7r 


k 2 dkP(k ) 


/ dy.A/"(y) / dxgf(k,x,x) 
'o Jo 


where the function g E is defined as 


9t(k,X,x) = 


XX 2 


(£ + 2 )\ 

(£- 2 )! 


T v (k,x) 


(6.43) 


(6.44) 


In Eq. (6.43), A/"(y) represents the distribution of sources as a function of the radial distance 
X defined such that Af(y)dy is the number of sources between x and y + dy. At order {0,1}, 
it is sufficient to consider the homogeneous source distribution, so that the observed shear 
and convergence for sources distributed up to y + are then defined by 


rx+ 


fX+ 


«o (x+,n 0 )= Af(x)^(x,n 0 )dx, 7oC X+,'ri 0 )= A^yb^y, n G )dy. (6.45) 


Since here J\f depends on y alone, this integration can be performed after the multipolar 
decomposition so that we perform the replacement, e.g. 


CX+ 


y)-> / A/"(y)E £m (y)dy , 


(6.46) 


in order to build the cosmological observables. Let us emphasize that this derivation can 
actually be performed in a simpler way [2]: since the source term derives from a potential, 
one could have simply used the Fourier transform directly in Eq. (6.41) and then expanded 


the exponential according to Eq. (5.9). The present derivation is however more general when 
used to higher orders {n,p}. 

On small angular scales, that is, in the limit l 1, it is possible to use the Limber 
approximation [83]. Such aproximation consists in using 


dxf(x)j e (x) ~ 



with L = £ + 1/2. If we commute the time integrals according to 



dy/(y,y) 




dy/(y, y) 


we arrive at the simple expression 


n EE{ 0 , 1 } 


l (l + 2)U 

4 (£-2)! Vl 


(6.47) 


(6.48) 


(6.49) 
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with 


2 


V, = f 


I o X" 


X 


rM 4,x 

X 


d%A/"(x) 


(x - x) 


XX 


(6.50) 


The angular power spectrum of the convergence k is obtained in a similar way. Indeed, 
if we consider only the dominant contribution of Eq. (6.42a) at small scales, it is sufficient 
to replace gf by 


9e( k ,X,x) = 


+ 1)X X h(kx)T*(k,x) 


(6.51) 


2 XX 

in the previous expressions to get . Using the Limber approximation, we then obtain 


^kk{ 0,1} 


4 


V, 


(6.52) 


and we check immediately that for large l, ~ cf E ^ 0,1 \ 

Finally, the angular power spectrum of the cross-correlations between the shear and the 
convergence is given by 


Cf { °’ 1} = IJ k 2 dkP(k) dxAf(x) d X9e(k,X,x)\ (J J d; X9e(k,x,X) 

for which the Limber approximation gives 


r nE{ 0,1} ^ £(£+ 1) I (i + 2)! 

1 ~ 4 


(6.53) 


E. Order {1,1} 

1. Geodesic equation 

In principle, we need to determine A; 0 ! 1 ’ 1 ) from the geodesic equation and then a;®! 1,1 !. 
As we shall see, these terms are only needed for the expression of the convergence 
We will instead focus on the computation of the cosmic shear and also the rotation 

VTM}, 

since they give the leading order of the 5-mode and the rotation. Fortunately, that 
computation does not require the solution of the geodesic equation up to order {1,1}. 


2. Sachs basis 

In order to get a definite expression involving only covariant and radial derivatives, we 
need to expand the direction vector n- around its background value n®, so as to use the 
, taking into account the contributions of order {0,1} and {1,0}, and 
similarly for the projection vectors n- a . We must use 

™-(x,<) = < + [cu a{0 ’ 1} (x,<) + w a{1 ’ 0} (x,<)] n° a ~, 

K(X, K) = rC - [w<°’ 1} (x, <) + ^’° } (x, <)] < ■ 


definition of 5 IV C 3 
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(6.54a) 

(6.54b) 












It turns out that only the expression for the projection vectors is needed since the direction 
vector n- appears only in terms which are already of order {1,1}. Additionally, we must 
convert the derivative along the tetrads ifj noted by di to derivatives along the Cartesian 
coordinates, and these are related from Eq. (4.8). This correction is only relevant for the 
term djpjcp because the other terms are already of order {1,1}. We thus use 


didjip = d.djip - 2 i B k ( l d j) d k (p , (6.55a) 

fa = I),D,B + 2 BSn + 2D(iBn°j + Bn°n °. (6.55b) 


3. General form 


Since in Eq. (4.34) the two terms *!o and ^1^7- Z ac V cb do not contain ’ 1 (because 

dk°/dx and 7 Z ab vanish at order {0, 0}), it can be integrated to give 


dV’S’ 11 _ g 

dx 


d In k° dV ab 1 

-i-; T - 9 7 Z ac T) cb 

dx dx (fc °) 2 


dx' • 


(6.56) 


(We remind the reader of our convention, in which we split k° and T> ab respectively as 
k° = —1 + 5k° and V ab = x$ab + d'D ab .) In the first term of the integral, given that dfc°/dx 
is at least of order {1,0} + {0,1}, the term D' ab can be expressed using the formulas found 
in the two previous sections, that is 


d ( yD'd’ t dA°- 1 } /"x /dfc 0, {i’0}/{0’i} 


dx 


dx 


-$ab + X7£ 


{1,0}/{0,1} 

ab 


dx- 


(6.57) 


Equation (6.56) can then be integrated as 


rad’bx) = 


X - X 0(1,1} 


X 


s ab’ (x)dx, 


where contains all source terms of order {1,1}. It is explicitely given by 


dh'(x) = 


\-R 

(. k °) 2 


ab „ xTi d ln /u° , , 0\ r , dfc° 

"I - XT^'acO^cb j _ (2 “h )0 a b “I - 


(6.58) 


d In x 


d In x 


dxx77ab + X 2 $x l di1Z ab , 


evaluated at order {1,1}, and where the last term arises from the fact that, at this order, 
there is a correction to be considered since we have to go beyond the Born approximation. 
That is, we cannot just integrate on the Friedmann-Lemaitre geodesic; instead we integrate 
on the geodesic x*(x, n 0 ) = X n o + &e*(Xj n o), so that the source term is 

Sab{x l {x, n o )) = Sab{x, K) + <fa*(x, n o)djS ab , (6.60) 

which implies that 

S { ab 1 } (x\x,n 0 )) = Sfi’ 1 } (xX) + fa 3 ’ { 1 - 0 } (x,n o )a 3 i°’ 1} , (6.61) 
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since 


d,S. 


{ 1 , 0 } 

ab 


0. It follows that the source term is explicitely given by 


+ 

+ X$ab 


n 


ab 


{1>1} d 5k°W . 

+ X -i- °ab 


(k°) 2 J ' x dy 

x («£%!'" + K<+ EEE £ *7& 0, d* 
dx 

+ x 2 te 3,{m djll { p i} . 


+ 


d(5fc 0{1,0} o,i } 


$X. 


xKr dx 


(6.62) 


We see on this expression that the general source S^’^Cx) has several contributions. First, 
it has contributions from the vector and tensor modes Bi and E %3 (respectively noted 
and S'^ 1,1 J' T ) which are at least of order {1,1} since they vanish in the pure FL case; they 
enter the terms 7and dk 0 ^ 1,1 ^. Then, all the other contributions are formally products of 
the scalar perturbations by the geometrical shear; they appear as products of {1, 0} x {0,1} 
terms. To compute explicitely these terms, we decompose the source term as 


cl 1 , 1 } 


(6.63) 


Each contribution can be further decomposed into its trace, symmetric traceless and anti¬ 
symmetric parts as 

Sab = 8abS + S{ab) + S[ab] . (6.64) 

Since our goal is to compute the effect of an anisotropic phase on the cosmic shear, and 
not on the convergence, we are mostly interested only in the symmetric traceless part. We 
shall thus not report the computation of the trace contribution to the trace part, except for 
the contribution coming from vectors and tensors, so as to be able to compare our results 
with the standard results in the literature, in the cases where the vector and tensor modes 
are considered even around a Friedmann-Lemaitre background. A full computation may be 
useful in order to cross-correlate weak-lensing with the magnitude of supernovae. 


4- Vector and tensor modes contributions 


The vector and tensor contributions are easily found from the literature [62( [83]. Splitting 
the vector mode into a radial and orthoradial parts as 


By Bi kjB r , Ejj Bij T 2f ^E^ T E r riVj , 


(6.65) 


the expression for the Riemann tensor given in Appendix B 2 for vector and tensor modes 
gives 


B { a i 1}T = nfn°£ \E”- - in q d [q E£ + n p n q (dyd 3 E pq + d p d q E.£ - 2n p n q d q d (t E j)p ](6M£ 
np 1}V = n£n£ + n q d q d {l B 3) - n q did 0 B q ) . (6.66b) 
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Using the projections of partial derivatives into radial and covariant derivatives (see Ap¬ 
pendix [ITT]), we deduce that the vector and tensor contributions to the sources (6.63) are 


5. 


{1,11V 

ab 


(x) = X5abB r , r + 


-D c D c B r + -^[ X 6 D c B { 


X 2 dx 


+ ~^(X 2 B r 
X dx 


+ — ( xD( a B b) ) - D {a D b) B r , 


(6.67a) 


s£’ 1}T (x) = xS, 


9E r 6 a b 

ab~i\ r 0 

dx 2 

d 2 

+ XX 7~2 (xB(ab)) H 


D c D c E r - X 2 ^~ 2 E r ~ 6 ^-( X E r ) - X 2 D C E C 

dx dx X dx 


dx 2 


D {a D b) E r 


2 d ^xE>( a E b) ), 


(6.67b) 


with the notation for the radial derivative >r = x l di. The hrst terms of each expression 
are respectively the V and T contribution of the term in dk 0 ^ 1,1 ^ in Eq. (6.62). For each 
of these two expressions, the hrst line contributes to the trace of the Jacobi matrix, that 
is, to convergence The second line contributes to the cosmic shear since it is 

symmetric and traceless. By construction there is no antisymmetric part, so the vectors and 
tensors do not contribute to the rotation . 

In order to compare and recover the results of Refs. [6 % H1J, we must use the fact that 
vector modes are transverse, and that tensor modes are transverse and traceless. This allows 
us to get the relations (see also appendix D 1) 


0 = D a B a + x(B r ), x + 2 B r , 
0 = D a E ab + 3E b + x(A'b),x ) 
0 = D a E a + 3 E r + x.{E r ), x ■ 


(6.68a) 

(6.68b) 

(6.68c) 


5. Trace free part of the quadratic contributions 


Starting from the general expression (6.62), the only contribution of the terms of order 
{1,1} to the trace free part is x 2 77yy Then, the terms are decomposed as 

= XK tm/mx b + n wwi (6.69) 

since, as we have just seen, there is no rotation at order {1,0} and {0,1}. To finish, it is 
obvious that 

d6k o{i,o}/{o,i} fdk°\ {m/{0 ’ 1} , , 

' . (6.70) 


dx 


dx/ 


It follows that the trace-free part of the quadratic contribution of the source term is given 
by 

{1,1} 


AA'Ax) = x 2 


n 


ab 


(k 


C02 


I v 2 7? { 1’ 0} k{°’1> 4- v 2 77 {0 ’ 1} k{ 1 ’°> 

+ If /C (a6> K + X ><-(ab) h 


(6.71) 


+X 2 K"«7 S’ + X 2 R {0 ' 1} 

I* 

+X 2 5x t{1 ' 0) d,K { °f . 




{ 1 , 0 } 
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All terms, except the first one, involve products of quantities which are {0,1} and {1, 0}, and 
have been already computed. Note that the first term is kept in the form (lZ/ ab \) / , 

and its detailed expression must be found using the perturbed Riemann tensor given in 
Appendix |B 2| This is indeed more convenient since we shall express everything in terms of 
the tetrad basis, and we will just need to use the fact that 


&r = A;-0 O U = (1 - $)A£. 


(6.72) 


We find for this first term 


77 


(ab) 


(k 


0}2 


{ 1 , 1 } 


i 3 
= 717 77 7\ 
(a b) 


@ij ffl -/■ 

n 


+ 2(p tr ) + 2 ak(fi-dj)ip - djpjtp - 
k 




77 


(6.73) 


+ 2 


( 


aikTl 

V 77 3 ) 77 ’ rr 


a u n k n l 0 <7i*n fc 

OiOj V + 2—r—dpt 


77 


u 


3 ^ ,r 


Then, we can split all partial derivatives into covariant and radial derivatives, using the 
expressions of Appendix D 1 This term is then given by 

{i,i} 


(lF)?) = -X 2 D { A } Z (2yv + V ') + 2 X D {a ED b)V 


(6.74) 


+2 X 2 D {a BD b) (?) + 2 D c D {a BD b) D c ip - (1 - 2 B)D {a D b) ip 

\ V / ,r 

9 /T r D (a D b )E\' 9 
- 2 a: 2 ( ’ Z ] +2y 2 


-A 




U 


U 




77 


A 


2-D(a-Dfc)£^ 


77 


-^Z7 (aJ D fc> T + ^D (a ED 6) ( - ) - 2y 2 D (a E' ro ^ 1} + 2y 2 ^ 


77 


77 




T> 


T {i,o} 

7 {a 


T>m ( ^ 


Finally, the last term of Eq. ( |6.71[ ) needs to be evaluated. It can be read directly from 
the previous results at order {1,0} and {0,1}. We need only to split it into radial and 
covariant derivatives using the formulas of Appendix D 1 We find that its contribution to 
the traceless part of the Jacobi matrix is given by 


= - 6 r^D {a D b) (v >r - 2^j 

-a c ^D c D {a D b) <p - 2 x 2 <*\l' 0 } D b) (H> 


(6.75) 


l (« 

Let us emphasize that, when a a ^ w a ) the source is partially seen on its side. 

To conclude, the source term (6.62) is obtained by combining the two terms (6.67) for the 
vector and tensor contribution to the {1,1} part, the term (6.73) for the quadratic scalar 
contribution and the term (6.75) for the post-Born approximation, to which we need to add 
the 6 terms which are products {1,0} x {0,1} in (6.71), obtained from the expressions of 
the former paragraphs. In principle, once all these contributions to the sources of the Sachs 
equation are identified and decomposed into radial and covariant derivatives, one should 
apply the formalism detailed in Appendix [Vj and expand each term in spherical harmonics 
for both the angular dependence and the Fourier dependence. 

This procedure is however extremely long and includes a large number of terms. We 
will not detail it here but instead just identify the dominant contribution and compute its 
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effect on the Jacobi matrix in order to derive the leading contribution to the 5-modes in 
the next Section. Indeed, once we convert the covariant derivatives into spin-raising and 
spin-lowering operators, each covariant derivative is clearly associated with a factor l. In the 
flat sky approximation, that is, in the small angle approximation, the dominant contribution 


arises from the first term on the second line of Eq. (6.75), 


h { ab 1} (x,n 0 ) ^ 


X-X-.cfi.oi 


Jo X 

as it enjoys three covariant derivatives. 


’ i (x,n 0 )D c D {a D b) (p(x,n 0 )dx 


(6.76) 


6. Trace part of the quadratic contributions 

As discussed in the previous paragraph, the computati on of the trace of involves 

a lot of terms such as T 0 ^ 1,1 ! and the fourth line of Eq. (6.62). This tedious computation 
can indeed be performed with all the details given in this article. It will however give only 
a small correction to k, the leading order of which is the standard convergence 

We thus decide not to include this computation here since we are mostly interested by 
the lowest order dominant effect related to the anisotropic expansion. 


7. Rotation quadratic contributions 


As we have seen, the rotation vanishes at orders {0,1} and {1,0} so that its leading 
order contribution appears at order {1,1}. Since 77^’^ i s symmetric, its contribution arises 
simply from the two first terms of the second line of Eq. (6.62), that is from the source term 


^H } (x) = X ( 


/C| a|c 0U c\b\ 


+ 7?{°’ 1 }A'n{ 1 ’°jA — 
T l^{ a \c 0U c\b\ 


)=e ab S^ 1} ( X )- ( 6 . 77 ) 


Using the expression of the previous sections, it is explicitly given by 

dUteO = -U+B+V r x(x z X) D.DMx)ix~ (~ ++ +) 

2 Jo X 

-'-D+D+tp [ X ( X-x)x L> _ D _ E / (~ )d ~ _(_<_>+). (6 .78a) 

2 Jo X 

The general expression for the rotation is then obtained through 

V {1 ’ 1} (x,n 0 )= f }( -^S { V o; l} {x, n Q )dx ■ ( 6 . 79 ) 

Jo X 

The rotation is thus sourced by the coupling between the usual cosmic shear of the standard 
scalar perturbation around a Friedmann-Lemaitre spacetime (D + D + ip) and the quadrupolar 
contribution due to the geometric shear (D^D^B" = TUTUS'). 
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8. Integration over the source distribution 


The last point that needs to be discussed before turning to the multipolar decomposition 
and the computation of the angular power spectra is the source distribution. 

The source distribution represents the mean number of object normalized to the mean 
density observed in a solid angle dfi 0 , that is 


dN 

dff 0 dy 




(6.80) 


In the Friedmann-Lemaitre and Bianchi I background spacetimes, which are both homoge¬ 
neous, J\f is constant on any constant time hypersurface, which means that it depends on y 
a lone Trims 

A/"(x, n Q ) = A/"(x) + A/” {0 ’ 1} (x, x l ) + A/" {1,0} (y) (6.81) 

where the second term is the standard fluctuation of the number density due to the large scale 
cosmological perturbations and for which it is understood that the position x l is evaluated 
on the background geodesic, that is x l = x n o- 

Note however than when one turns to redshift space, on which the observations are 
actually performed, one needs to take into account that z is a function of y and the direction 
of observation n Q , so that we should rather use 


A f{z,n Q ) 


= a r( z ) + a r {0 ' 1} { z , x l ) + a r {m (z, n Q ) , 


(6.82) 


where again the position is evaluated on the background geodesic with x l = x n l- If follows 
that, when computing the observed quantities, 


l { ab l} {x,n 0 ) = / dy [AT(y) 7 I fc 1 ’ 1} +A/' {0 ’ 1} (y,a: i )7 ( [ b 1 ’ 0} +AA 1 ’°I(y) 7 I fc °’ 1} 


The second term is the standard correlation between the fluctuations of the source distribu¬ 
tion and the cosmic shear. It inherits a directional dependence from the spatial dependence 
A/’I° ,1 I(y, x l ) given that it is evaluated on the background line of sight, that is with x l = x n o- 
Because of the coupling to the pure ZTmode 7 ^’°^ it will induce IT modes in the source av¬ 
eraged cosmic shear. This component is expected to be important on large angular scales. 
The third term is a correction that arises from the fact that the formation of structure differs 
a priori in the presence of a geometrical shear, but it does not contribute the H-modes since 
it does not have a directional dependence. However, it induces a correction for the IT modes 
and for the convergence. 

Now, in redshift space, one needs to be more careful since 


{l.i} 
Tab T 


z, n, 


= / dz 


A/"(5) 7 2’ 1} +Af {0 ’ 1} (z,x i )jH fi} +A/' {1 ’ 0} (S,n o ) 7 2’ 1} 


Both the second and the third term depend now explicitely on the direction of observation, so 
that the convolution by the source distribution has to be performed before the decomposition 
in spherical harmonics, and both terms will generate H-modes out of the ITniodes of 7 ^’°^ 
and respectively. However, these effects should not dominate for small angular scales 

and we shall thus neglect them. 
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VII. ORDERS OF MAGNITUDE 


The previous sections have provided all the elements needed to compute the contribution 
of the 5-modes at order {1,1} and their correlations with the 5-modes and the cosmic 
shear. It is obvious that any further computation has to be performed numerically. It is 
however important to exhibit the dominant contribution. 


A. Dominant effects 


Once the covariant derivatives are expressed in terms of spin-raising and spin-lowering 
operators, it is rather straightforward to realize that any covariant derivative is associated 
with a factor £ in multipole space. The dominant terms contributing to the shear are thus 
the ones with the highest number of covariant derivatives applied to p. 


For instance, at order {0,1}, the convergence is dominated by the last term of Eq. (6.42a) 
on small scales. That is 


hi 


{o,i} 1(1 + 1) f x X-X 


lm 


XX 


<Ptm(x) dx, 


simply because of the geometrical factor £ 2 . It is indeed the term which is usually presented 
in textbooks. This term dominates over the second one even at small £, i.e. for £ > 2 — 3, 
that is for all practical purposes. 

When applying this small scale approximation scheme at order {1,1}, we realize that 
there is just one dominant term - the first one on the second line of Eq. ( |6.75 ) - which 
possesses three covariant derivatives, that is 


h 2’ 1} (X', n Q ) ~ - f ^-W^a c{1 ’ 0} (x, n 0 )D c D {a D b) ip(x, n Q )dy . 

Jo X 

Physically it corresponds to the orthoradial displacement of the central geodesic on which 
the Sachs equation is evaluated, when compared with the background geodesic. It is as if the 
sources of order {0,1} contributing to the Jacobi map had been lcnsed by the orthoradial 
displacement of order {1, 0}, resulting in an order {1,1} effect. This is similar to the lensing 
of first order sources of CMB around the last-scattering surface by first order gravitational 
potential in the foreground, resulting in a second order lensing effect in the CMB. 

The first consequence of this is that the formalism used to compute the CMB lensing can 
also be applied to obtain the resulting Jacobi map due to this leading order term. However, 
there is a slight difference. Indeed, for the CMB the sources are all located in a background 
around the last-scattering surface, for which there is a deflection due to the gravitational 
potential crossed between emission and reception. For the general solution giving the Jacobi 
map, however, the sources are distributed from the observer up to the maximum redshift 
of the survey. For each source there is a different deflection angle as it depends on the 
trajectory between the source and the observer. 

Finally, we must remind that the treatment of CMB lensing by a gradient expansion [85] 
holds only until the deflection angle is comparable to the angle of structures in the CMB. 
Beyond that scale, this method underestimates the effect of lensing and one has to resort 
to a full-lcnsing method as exposed in Ref. [7SJ [86j. Since we are interested in an order of 
magnitude estimate of the effect of geometrical shear on the cosmic shear, we will present in 
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the next section a gradient expansion method based on Ref. [85], but one should be aware 
that for any amplitude of the geometric shear, there must exist a scale £* beyond which 
this treatment is inaccurate. The method for the full-lensing method is exposed briefly in 
Appendix [Cj 


B. Lensing of the central geodesic 


1. General formalism of the gradient expansion 

Any observable at a given affine parameter y hi a given direction n Q is formally obtained 


from an integration on the background geodesic over its sources given by Eq. (5.4), that is 

X*(x,n 0 )m s Q = f S x (x,X,n 0 )m s 0 dy. (7.1) 

Jo 

However, and as discussed above, a true observable like the cosmic shear is obtained by 
averaging over the true normalized profile A7(y) of sources as 

poo poo px 

X s (n 0 )m s 0 = dyA7(y)X s (y,n 0 )m*dy= / A7(y)dy / S x (x,X,n 0 )m a 0 dy. (7.2) 

Jo Jo Jo 

Note that the integrals can be interchanged using 

poo px /»oo /»oo 

/ dyA f(x) d xf(x,x)= dy / dy A7(y)/(y, y). 

Jo Jo Jo J\ 


(7.3) 


We consider only the effect of the dominant term in Eq.(6.75), which corresponds to the 


lensing of the sources, that is, it transforms the sources according to a parallel transport 
along the lensing vector «. A lensed observable A" s is then obtained from an integration 
over the lensed sources. If the lensing effect is small, it is sufficient to use a Taylor expansion 
of the lensed sources to express them in terms of the unlensed sources, the small parameter 
being the lensing vector ol. Furthermore, if the lensed vector can be written as the gradient 
of a scalar, as a a = D a a , then at lowest order in the Taylor expansion, we get for the lensed 
source 

S x " (x.x.n°)"*s = S x '(x,X,n 0 )m’ o + D a a(x,n<,)D a [S x ‘(x,X, n„)m‘] . (7.4) 

(7.5) 


Using Eq. (D26), the multipoles are easily extracted as 

~A _ oX s 


S£(X,*) = S£(X,*)+ E , 

7711 ,^ 2,7712 

where the are the multipoles of the lensing scalar when decomposed into spherical 


harmonics, and the coefficients m ' 2 are defined in Eq. (D28) 


2. Multipoles of the lensing vector 

The previous expression depends on the multipoles of the lensing scalar, that can actually 
be obtained very easily. First, following the definitions (6.5 6.6) we define a matrix such 
that 

a±(w 0 , X) = D±a(n 0 , y) = D± -a^ y)n*< 


(7.6) 
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Given Eqs. (6.13b) and (6.75), the components of ap (x) are just 




-/%( 0) + 2 [ dx'~ — —o-ij(x') 
Jo X 

Pij (°) - - f dx'Pijix ') ■ 

X Jo 


(7.7a) 

(7.7b) 


Then, similarly to the computation of the coefficients £ 2 m in Eq. (4.26), the multipoles of 
a(n 0 ,x), defined by a = X) m a 2 mb 2m reduce to a quadrupole and are explicitly given by 


«2o(x) 



(x) + a yy (x)\, 


Q; 2±2(X) 



a yy(x)] ) 


(7.8) 


if the coordinates system is adapted to the eigendirections of the geometrical shear. 


3. Extracting the spatial shear components from off-diagonal correlations 


A byproduct of the formalism just developed is that we can extract information about 
the geometric shear ay,- from cross-correlations between the E- and 5-modes multipoles of 
the cosmic shear, E^ m and 5^ m , and the multipoles of the convergence n, that would 
otherwise vanish in the pure Friedmann-Lemaitre case. Indeed, even if the 5-modes are not 
sourced initially, as is the case of a Friedmann-Lemaitre background, at the perturbative 
level there will be a lensed 5-mode term sourced by the 5-modes of the background shear. 
In order to extract this effect we decompose the (lensed) 5- and 5-modes of the source as 


s tn (x, x) = Sf m (x, X) ± x) 


(7.9) 


with a similar decomposition for the (unlensed) Sf m . 
follows that 


Then, using the properties (D28), it 


s?M,x) = -i £ «2 mt (x)sZ m (x,x) 


mm pm2 


mi,m2 
£ 2 =£±1 


sFJx,x) = sfjx,x)+ £ «2 TO K)s£ 2 U,U + 2/Sr™. 


(7.10a) 

(7.10b) 


mi,m2 

£ 2 = e /±2 


We remind that there is no tilde on Sf m on the right-hand side of the above equation, since 
it corresponds to the unlensed sources. Since the convergence is a spin 0 quantity, then from 


Eq. (7.5), its sources are transformed under lensing as 


SL(X, x) = sux,x)+ £ c* 2mt (x)s t ~ ma (x, x) G” 


mm\rri2 


(7.11) 


mi,m 2 

e 2 =e r £±2 


From these expressions, it is clear that the off-diagonal terms coming from the 55, 55, kk, 
k5, nB cross-correlation matrices allow us to put constraints on a 2m and, consequently, on 


the geometric shear components by means of Eqs. (7.7). To see how that is possible, we 


must remember that the sources should be convolved with M(x) by means of Eq. (7.2). 
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In Section VID 5 the EE, kk, and kE correlations at order {0,1} (that is, without the 
effect of tensing by the geometric shear) have been computed and they are of the form 


Cf z = 


dxi / dxi7V(x!) / d\2 / dy 2 7V(y 2 )C'f (Xi,X2,Xi,X2) (7.12) 


'Xl 


'X2 


where the indices X and Z take the values k, E, and with the source correlations 

(Xl j Xl ')^£ 2 m2 (^2; X2)} ^t\£ 2 ^m\m 2 ^£^ (Xl j X2j Xl i X2) j ( 7 . 13 a) 


Cf Z (Xi,X2,Xi,X2) = - 

^ Jo 

For the lensed observables, we define similarly 


k 2 dkP(k)g^ (k, Xi, Xi)9e (&, X 2 , X 2 ) • (7.13b) 


' £M ~ 


dxi / dxi7V(yi) / dy 2 / dy 2 A7 (x 2 )« 2 m(Xi)Q(X i, X 2 , Xi, X 2 ) 


'Xl 


' X2 


2 

7T 


fc 2 dfcP(fc) 


dyi / dx 1 A7(y 1 )a 2M (x 1 )^f (fc,xi,xi) 


x 


0 

dy 2 


'Xl 


dy 2 A/" (x 2 )gf (k, x 2 , X 2 ) 


' X2 


such that the following non-vanishing correlations are expressed as 

/R 7?* \ i jmM(m—M) ^yEE 

\ D im^t±\m-MI ~ 1 +‘2 1 £2e±l 'tilJlfi 

/ r \ ; jmM (m—M) -jyEn 

\ n tm^£±l m -Ml ~ 1 + 2 V 2 fcti 


(7.14a) 


(7.15a) 

(7.15b) 


Not only do we get off-diagonal contributions for P-modes with the P-modes and the con¬ 
vergence, but we also get off-diagonal correlations between k and P-modes. They read 


{EcmEi-I2m- M ) — +^£2£±2 r l±2M 


T mM (m-M)spEE 

• I 


\M 


(n£m K i±2m-M) ~ ^£2£±2 
\^£m K £±2m-M) = +‘2 1 £ 2 £± 2 


m.M(m-M)~j KK , -i\Mt 

>i±2 m + -U h 


* \ _ jmM (m-M) 'ryEn 

1 0 o /?_lo r , 


+ ( 1) +2 1 e±2 2e 

7 i—M 
£±22 £ 

£±2M + ( — 1) ± £±22i 


j(m—M) —M m^yEE 
7 0 -lo off ' £ M 1 


M T (m—M) -Mm^KK 

' £M 1 


1 /_i\M -rynE 

“r t ^ J 7 ^4-0 7 c / 1 


IM ■ 


(7.16a) 

(7.16b) 

(7.16c) 


Note that in all these cross-correlators, M ranges from —2 to 2, thus spanning the five 
degrees of freedom of the lensing potential a 2 M, and consequently of the underlying Bianchi 
geometrical shear ay,-. These expressions for the correlators are however not ideal to relate 
the correlations to the lensing potential, and thus to the components of ay. Instead, we 
define appropriate combinations of the correlators by resumming them as 


- E X-D 


m+£ 1 +£ 2 


Ci 2 t 2 
-m M m — M 


) (X,.„z ( * , m _ M >. ( 7 . 17 ) 


For instance, for the EB and EE correlations, we get 
BE j^2M — ST' , /S7_1 )m+l 


H£± 1 


EE a2M _ 


l£± 2 


m 


e 2 e±i . 

-m M m-M ) \^™^±i , m - M ) > 


«■ 2 l±2 * > 

-171 M m - M J • 


(7.18a) 

(7.18b) 
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Then, by using the definition of the symbols + 2 ^i ±™' and I^KT' an d the orthogonality 
relations of the Wigner 3j symbols, we get 


'A 

,J A 

* A 

* A- 


2 M 
W± 1 

■ 2Fm±l ^EE 

- 1 y /5 e±1M ’ 

(7.19a) 

2 M 
«±1 

;Fm±l r,E k 

— 1 yjT r l±\M i 

(7.19b) 

2 M 
M±2 

2Fm±2 y) EE , 2F(.±22i y) EE 

~ 1 L±2M ' ' IM J 

(7.19c) 

2 M 
v ll±2 

Fm±2^, KK , F ( .±22i /T)KK 

— ' l±2M "r ^ r im , 

(7.19d) 

2 M 
l U±2 

2Fg2£±2 ^Ek , F t ±22i ^ kE 

~ ^ ' t±2M ' r IM 1 

(7.19e) 


where the symbols 2 Fn x t 2 are defined in Appendix |D 4[ 

Approximate expressions of this correlators can be obtained in the Limber approxima¬ 
tion (6.47) and read 


’£ M — - ~ A - VlM, 


t> kE ~ — 

'PM — 'PM — 


Ek, 


IM 


IM 


1(1 + 1 ) 1(1 + 2 )! 
(( - 2 )! 


V. 


‘ M 


~ 1 (1 + 2 )! 

,M ~ 4 (1-2)1 ' 


IM 


(7.20a) 

(7.20b) 

(7.20c) 


with the function Vim given by 


Vim = 


dx 


'0 X 


X 


~^P ~ «2 M(X) 


TV[ ~A 


dyAT(x) 


(x - x) 


XX 


(7.21) 


where we used the notation L = £ + 1/2. 

This provides the expressions of the 5 (off-diagonal) correlators (7.19), each having 5 
components, and all being linear in o iy We stress that the measurement of these quantities 
from further surveys will allow us to get stronger constraints on the spatial isotropy of the 
universe, thus pushing forward the “beyond ACDM” program. 


4■ Autocorrelations of B-modes from the lensing of the central geodesic 


The previous off-diagonal correlators are the most direct consequence of a late time 
geometrical shear on weak-lensing. However, experiments are mostly designed to measure 
the diagonal part. In this section we compute the autocorrelation of H-modes induced by 
the dominant lensing term. This angular power spectrum will thus be quadratic in cr^. 
Contrary to the previous estimators, it does not allow us to reconstruct the full geometrical 
shear a t] , but can be used to set constraints on a 2 . 


Using the properties of the Wigner 3j symbols given in Appendix D 4 and starting from 
the lens sources (7.10), we obtain that the H-modes angular power spectrum of weak-lensing 
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cosmic shear generated by the lensing of the central geodesic is 


r°° r x 1 r X2 

cf {1,1} = / dxi / dx2ftf(xiW(X2) / dx'i / dyv 


'0 


x 


Ct 2m(X / l)«2m(X2) r yE£;^ / / n, 

2^ -c-Q+s(Xl,X2,Xl,X 2 ) 


/ , s\2^m+s 


Fp 9/4-.e ) 2 


s=d=l,m 

If we now factorize the time integrals we simply get 


2£ + 1 


(7.22) 


c bb { m } _ 2 r° 
_ 5tt/o 


fc 2 dfcP(fc) ^ 


s=±l 


f2-Pl2l+sJ 

2 £+ 1 


E 


Note that for large £, the P-coefficients behave as 


dxN(x) / dx , a2m(x')ge+s(k,x,x' 

Jo 

(7.23) 


Fp 9 /-i_i) 2 15 ( 2 -P? 2 i ) 2 15 


U- r ^2£+ly — ,. 

lim --- = — . hm 


£-> oo 2£ + 1 


2t r 


£—>oo 2£ + 1 


2tt 


(7.24) 


Apart from the six terms (sum over m = —2,0,+2 and over s = ±1), this is numerically 
as fast as computing the correlation Cf E at order {0,1}. Then, using the Limber approx¬ 
imation (6.47), with the definitions L s = L + s = £ + 1/2 + s and £ s = £ + s, it leads 
to 


a 


BB{ 1,1} 


1 r d X Y- {is + 2)! 0 / L„\ ( 2 Fp 2 p +s ) 


20 


E 


V f-, X - 2 )\ V X J 2 f+l 


P ~ 


X 


E 


X 


T* ( - 2 ,x I «2 mix) / dxM{x) 


(X ~ x) 


XX 


(7.25) 


Finally, using (7.8) to get 


V |a 2m (x)| 2 = j^OdjixW^x) = j^Wx )\ 2 

m 

we obtain a very compact expression, valid only for large £, which is 


(7.26) 


c bb { i,i } *x P (i 

1 10 J 0 X 2 \X 


X 


T* l-,x «(X) / d X A fix) 


(x - x) 


XX 


(7.27) 


A numerical analysis of some simple anisotropic phenomenological models, together with 
observational constraints from Euclid [7l] and SKA [?2j surveys, will appear in a companion 
paper [70] , 


VIII. DISCUSSION 

In this article we have derived the observational signature of a late time anisotropic ex¬ 
pansion on the weak-lensing observables. To that purpose, we have provided all the technical 
tools, including the evolution of the background spacetime, the perturbation theory, the de¬ 
scription of the evolution of a geodesic bundle and the manipulation of observables on the 
celestial sphere. 
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Our strategy is to adopt an approach based on the observer point of view, in which all 
observables are expressed in terms of the direction of observation at the observer. Since 
a full solution to the problem cannot be attained straightforwardly, and given that CMB 
observations suggest that spatial anisotropy cannot be too large, we have developed a small 
shear approximation scheme. It allowed us to identify the following contributions compared 
to the standard Friedmann-Lemaitre case: 

1. the tensor and vector contributions to the source of the Sachs equation, which starts 
at order {1,1}, and the contribution of the scalar modes coupled to the geometrical 
shear, which is of order {1,1} as well; 

2. the evolution of all the perturbative modes, that is, of the transfer functions, which 
are decomposed as T Xs (k,t), where the dependence with the direction of k comes 
from the coupling with the geometrical shear in the Einstein equation (for that, see 
the Appendix [A]) ; 

3. the fact that the geodesic deviates from is Friedmannian form and which leads to 
post-Born corrections; 

4. the effect of the source distributions which are affected by the background shear or 
the scalar perturbations - that are respectively at orders {1,0} and {0,1} - and for 
which we would in principle need a theory of structure formation; 

We have then argued that the dominant term is related to the orthoradial displacement 
of the central geodesic on which the Sachs equation is evaluated, when compared with the 
background geodesic. 

While we have provided all the elements to perform the full computation, we have focused 
on this dominant term and demonstrated that there exist 5 off-diagonal correlators between 
Ei m , B em and K£ m each of which has 5 independent components and thus allow in principle 
to fully reconstruct the geometrical shear a t j. All of them are linear in a/7~L and only two of 
them involve the 5-modes. We advocate that their measurements in future surveys such as 
Euclid and SKA, on scales where the linear regime holds, can set strong constraints on the 
anisotropy. The amplitude for these two surveys is estimated in our companion article m- 

The existence of non-vanishing 5-modes also reflects iteslf in the existence of an angular 
power spectrum that is quadratic in <7/5. While probably easier to measure, it does not 
allow one to fully reconstruct the shear cry. 

This analysis sets the ground for stronger constraints on an anisotropic expansion, and 
possibly on the anisotropic stress on the dark energy sector. The new estimators that we 
proposed will also allow the control of systematics, and are new in the weak-lensing literature. 
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Appendix A: Perturbation theory in Bianchi / universes 


This section summarizes the general framework of linear perturbation theory in Bianchi / 
universes. Our approach is an extension of the formalism we introduced in Ref. [22], where 
perturbation theory was developed in the context of an early anisotropic stage. Here, we 
adapt this formalism for the physics of the late-time universe. Before we introduce the 
parameterization of the perturbations and the whole machinery of gauge-ivariant linear 
perturbation theory, we summarize some results regarding the appropriate Fourier transform 
in anisotropic spacetimes, and show how they can be used to extract the dynamics of scalar, 
vector and tensor modes from Einstein’s equations. We then use these tools to decompose 
the background shear and anisotropic stress in a general basis adapted to our coordinate 
system. 


1. Mode decomposition 

a. Fourier and SVT decomposition 

In order to correctly describe the dynamics of perturbative modes one needs a complete 
set of spatial eigenfunctions adapted to the symmetries of the spacetime one is dealling 
with. Since Bianchi / universes are spatially flat, at each constant time hypersurface these 
eigenfunctions are standard plane waves. Therefore, any scalar function of the comoving 
coordinates {x 1 } and time can be Fourier decomposed as 

/ a 3 u. „ 

J^f(ki,v)e ikiX '■ (Al) 

with the obvious inverse transfomation. Due to the lack of rotational symmetry, the direction 
of a wave vector will vary with time. In particular, since ki is constant, k l = rf^kj vary with 
time - its rate of change being given by: 

(% = -2cr% . (A2) 

Note however that kiX 1 = k l Xi remains constant. From the above expression, one can easily 
deduce the time evolution of the modulus k 2 = &%• and unit vector k l = k l jk as 

j = -a ij kikj , (fry = (cr jl k j ki)k i - 2(7% . (A3) 

As we are going to see, these expressions are crucial for extracting different perturbative 
modes from Einstein equations. 

Once equiped with a Fourier transform, we can proceed and decompose any 3-dimensional 
geometrical object in terms of its scalar, vector and tensor pieces. We start by decomposing 
any (3-dimensional) vector V) in its longitudinal and transverse pieces as V) = d t V + V^, with 
d l V % = 0. In Fourier space, this decomposition is equivalent to 2 

V i = k i V + V i , k% = 0. (A4) 


2 Note that we can always reabsorb i factors in the terms of the decomposition. 


59 




Since V % is orthogonal to k l , it can be further decomposed as 


Vi=Y, V a(k\v)e a (ki) 
0 = 1,2 


where {e“} represents a 2-dimensional basis defined so that 



= S' 


ab 


(AS) 


(A6) 


Analogously, any (3-dimensional) symmetric tensor can be decomposed into a trace 
plus a traceless part as V l3 = Tj i3 + A i3 S + 2 d^Vj) + 2Vjj, where A ?;j = d,dj — r y i jA/3, V) is 
transverse and Vjj is transverse and traceless. In Fourier space, such decomposition becomes 


Vij — T^/ij + ( kikj —^ ] S + 2 k(jVj) + 


(AT) 


where V) is given by Eq. (A5). V l3 is a transverse and traceless tensor decomposed as 


Vi]= Y. V A‘,I 

A=+, x 


(AS) 


with the traceless (e^-y*- 7 = 0), transverse (e^k 1 = 0) and perpendicular (e^-e^ = 5*) polar¬ 
ization tensor being defined as: 



11 _ 2 2 12 , 2 1 
i j 1 j -5l+ * j V i j 5* 


V2 


y/2 


(A9) 


Given the above decomposition, the correspondence between SVT components of any 
geometrical equation can be extracted uniquelly. For example, the scalar part of any vectorial 
equation of the form V) = 0 can be extracted by projecting it along k l , whereas its vector 
part can be extracted with the help of the projector 


Pij = la ~ hkj = e\e) + '7''). 


(A10) 


Likewise, the scalar components of any tensorial equation like V tJ = 0 can be extracted by 
projecting it along y^- and T t j , with the later projector defined as 


Tij — k^kj g7y 


(All) 


The remaining vector and tensor pieces can be extracted with the help of PfkJ and A“j\ 
respectively, where 

A ab = papb _ \p. jP ab ( Al2 ) 

In conclusion, the SVT degrees of freedom of any tri-dimensional vector and tensor are given 
explicitly by 


Vij 


k j Vj ) ki + PfVj , 


Ay kl Vki ) lii + ( i T« l V kl ) T,j + 2fc(. 


kh 


P m k n V 

r j) ^ v rr 


4 - \ mn v 

1 1 v mn • 
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b. Mode evolution 


The SVT decomposition introduced above is based on the properties of a tensor according 
to the rotation symmetries of the background spacetime. As such, in Friedmann-Lemaitre 
spacetimes this decomposition is always possible and will hold during the entire cosmic 
evolution. In the Bianchi / case, on the other hand, this decomposition will hold, strictly 
speaking, only on a given constant-time hypersurface. Due to the anisotropic evolution 
of space, SVT modes which are initially decoupled will couple nontrivialy as time evolves, 
implying in a set of coupled dynamical equations already at linear order in perturbations. 
Therefore, it is important to have expressions for the time evolution of basis vectors and 
polarization tensors, which will be directly dep endent on the spacetime shear. We have 
already met the time evolution of k l ; Eq. (A3). For completeness, we also give the time 
evolution of the vector ef 


and polarization tensor e(- 


(k 1 )' = (a jl k j k l )k i ~2a ij k j , 
( e aY = ~ > 


(4) 7 = - 4 M 4) p ij - (o*Pk) 4 + 4 44a 


(Al4a) 

(Al4b) 

(Al4c) 


Special care is needed when extracting SVT modes from Einstein equations, for the 
projections of SVT modes do not commute with time evolution anymore. As an illustration, 
let us consider the extraction of the scalar component of an equation like (Vi)' + "HVj = 0, 
where V) is any transverse tensor. In FLR W this equation does not have a scalar component, 
since k l Vi = 0. However, due to Eq. (A3) we now get 


k‘ [(V 4 )' + W] = (k‘V ,,)' - (k‘yV, = 2a%V, 


(A15) 


which is only zero when o tJ = 0. Further mode-extracting relations can be easily found in 
an analogous manner. For a comprehensive list of relations the reader can check Ref. 


c. Background shear and anisotropic stress 

Both the (background) spatial shear cr ?J - and spatial anisotropic stress n ?J are transverse 
and traceless tensors. As such, each of them is described by five independent degrees of 
freedom, which are best described in the basis {fci,e|,ef} adapted to the modes we are 
considering. In this basis, these two tensors can be written in terms of 10 new scalar 
functions as 

Vij = a n + 2 °V« 4 e “) + °ta 4 ’ (Al6a) 

np = — (kjtj — —7 ij \ Hu + 2 4^ n v a k(i e j) + 4^ n x A . (Al6b) 

A ' a=l,2 A=+,x 

It is important to note that these new functions, which depend of both ki and time, are not 
the Fourier transform of their respective tensors, which in fact are homogeneous and depend 
only on time. In other words, the dependence of (cr (| , a va , cr T \) and (n ., n va , H t a) with ki 
arises solely from the local anisotropy of space. 
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With the help of background Einstein equations (2.18a) and (2.18b) written in conformal 
time and the mode evolution described by Eqs. (A14a), one can show with a bit of algebra 
that 


(On = a 'u + ma \\ + 2 X77 = Ka2n n » 


(Al7a) 


v' _ 

)y a = (T ya A 2'B(Jya ~cy^\ a ^\\ A ^ ^ aAI^ ; 


b, A 


(ct')tA = v' T \ + 27/^a - 2 ^ M^ b a v a<J yb = kq 2 fl TA , 


(A17b) 


(A17c) 


a,b 


where A 4* b i s defined as [ 22 J 


M 


A 

ab 


l 

71 


1 0 
0 -1 




2. Gauge-invariant variables 


a. Geometry 

The most general linearly perturbed metric over a Bianchi / spacetime can be parame¬ 
terized as follows 


ds 2 = a 2 [— (1 + 2 A)dr ^ 2 + 2Bidx l dr} + ( 7 \j + hjj)dx l dx 3 ] (A18) 

where A is a free scalar function and 


Bi = diB + Bi, (Al9a) 

hij = 2 C ^ 7 ij + A 2 didjE + 2d^Ej) + 2 _ZT^-, (Al9b) 

defined together with the usual transversality and trace-free conditions: 

diB* = 0 = d t E\ E\ = 0 = diE ij . (A20) 


Under an active coordinate transformation, the coordinates of any point will change 
according to 


X*-> S'* = X* - £**(: x v ) , 

where the gauge vector is itself decomposed as 

(A21) 

e = T, e = d l L + E 

(A22) 


with diL 1 = 0. Under the transformation of Eq. (A21), the perturbations of the metric will 
transform as 


Sg^iu —> 


(A23) 
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where C^g^v is the Lie derivative of the background metric along £. Using the above param¬ 
eterization and the mode decomposition introduced in § A1, it is straightforward to show 


that the scalar and vector metric potentials transform respectively as 


A -7 A + T' + HT, 

(A24a) 

B^B T + {k v , 

(A24b) 

C -7 C + HT , 

(A24c) 

E —y E + L, 

(A24d) 

Bi -7 /I, - ',p/,') , -2i/Urr 0 />'/.. 

(A24e) 

Ei —y Ei + Li , 

(A24f) 


whereas E %3 is automatically gauge invariant. Based on these transformations, we can con¬ 
struct the following gauge-invariant variables: 


$ 

vL 


A -f- - s a 
a 

-c-n 


B 


B 


(. k 2 EY 
" k 2 
(. k 2 Ey 
k 2 


B,-~ :ij (E i )'~Ak i cTi i I>’E. 


(A25a) 

(A25b) 

(A25c) 


It is easily verifiable that whenever a = 0 the Fourier wave vector k will be constant and the 
above variables become the standard Bardeen variables for a Friedmann-Lemaitre universe. 


b. Matter sector 


Moving forward, we now parameterize the perturbations of the energy-momentum tensor 

(A26) 


defined in Eq. (2.14), which can be decomposed as 

T,f = + ST? 


From the normalization condition of the fluid total four velocity we can write 


Su M = -(-A, v ®), v l = d l v + v l . 
a 


(A27) 


with diV 1 = 0, as usual. Likewise, the perturbations to the energy density (Sp), pressure 
(SP) and anisotropic stress (Stt)) are introduced as follows: 


5T 0 ° = -Sp, 

5T° = [p(l -«-)^_lb,;(vd-/^). 
5Tq = -p( 1 + w)v l + Y J 57r j0 , 

ST) = 5PS) + 7 a «57ry - n jk h 


ki 


(A28a) 

(A28b) 

(A28c) 

(A28d) 


where Bi and h % j were defined in Eqs. (A19). Special care to the notation is in order here 
because, as one can check, g^ x S7 ^ Sn% = 9^ X ^\v + bl A v Sg x>l . 
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We also need to parametrize the perturbed anisotropic stress tensor Sir*. From the 
transversality condition ( u M + 5u tl )(Yl^ u + 8n lw ) = 0, we conclude that 

5tt 00 = 0 , 5n 0i = -II ^v 3 . (A29) 

Note however that these conditions do not fix We therefore further decompose Sn t j as 

Snij = 2 [n T 'y ij + didj7i s + <%7rj^ + 7r^] (A30) 


where T in ir T stands for “trace” and where, as usual, we have 


avf = o = <9v; 


•ijr 
ij ’ 


Ti n 
7T „• = (J . 


(A31) 


Moreover, note that ^ 0, which is why the above decomposition tensor 

has a trace. 

Under the gauge transformation (A21) and using again the appropriate Fourier decom¬ 
position, the above variables transform as 

5p —> 5p + p'T , 

6P -> 5P + P'T , 

(ik 2 Ly 


v —> v — 


v -)• v 


k 2 ’ 

(V)’ + 2i k j aijP li L . 


(A32a) 

(A32b) 

(A32c) 

(A32d) 


These transformations suggest the introduction of the following gauge-invariant variables 


Sp = 5p + p' 
<yp = 5p + p' 


B 


B- 


(k 2 E)' 

k 2 

(fc 2 uy 

k 2 


v 

v l 


= v + 


(k 2 E)' 


k 2 
E +BP 


(A33a) 

(A33b) 

(A33c) 

(A33d) 


The perturbed variables in Eq. (A30), on the other hand, do not have simple transforma¬ 
tions as above, essentially because there is no simplifying relation between the background 
tensor fl^ and the wave vector k l . Using 


An# = nhr + 

we find that 

^ -► w T +(-\T% + ^ i m i ^T, 

77 s ^ 7T S - ^r 2 T lJ % T + n ijtiktL - m, 

-)• ttY - T + ikPiUj&L + P’W-,1 . 1 , 

_T , T | ^ a^TT' T 1 

W ij + 2 P lrn 1 > 


(A34a) 

(A34b) 

(A34c) 

(A34d) 
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where Pij, T\j and were defined in Eqs. (A10), (All) and (A12). From the variables 
above we construct the following new variables 


7T T = TT T + 


- ''/'Ml', + la ij H 

4 2 3 


7T S = 7T S - '/'Ml',. B 


4 k 2 l] 


V 
2 jp\l 


B- 


(k 2 E) 

k 2 


7T, 


y = vrf - ^P?k l n' ( 5 - 


k 1 jl 


*1 = 7 xl + ^Arni , I B - 


IrriT t / 


U 




(k 2 E) 
~E 2 
[k 2 E)'\ 

) 

{,k 2 E)' 
k 2 


k l 


— UijEk 3 E + iIT,j — F- 7 , 

K> 

ikpm.UE - /yii, 


(A35a) 

(A35b) 

(A35c) 

(A35d) 


which, as one can easily check, are gauge-invariant. 


c. Gauge choice 

From the construction of gauge-invariant variables presented above, it is clear that an 
enormous simplification will be achieved if we work in a gauge where 

B = E = 0 = E i . (A36) 

In this gauge the scalar modes become 

$ = A , ^ ——C, Sp = 5p, 5P = 5P , v — v, tt t = n T , n s = n s , (A37) 

whereas the vector and tensor variables become 

= Vi = Vi + Bi , = 71 -Y , 7T5 = 7rg . (A38) 

Apart from the spatial velocity v, tl in this gauge the gauge-invariant variables coincide with 
the original potentials. In other words, by working in this gauge the final equations can be 
trivially (again, apart from v t ) replaced with the same equations satisfied by gauge-invariant 
variables. Note that this choice fixes the gauge completely, and is slightly different from the 
choice made in [ 22] , 


3. Perturbed Einstein’s equations in Bianchi I 


We have now everything needed to obtain the fully mode-projected and gauge-invariant 
Einstein equations. This is a tedious but straightforward procedure which requires careful 
comutation of time derivatives and the Fourier vectors through the use of Eqs. (Al4a Al4c). 
We note that the main difference with the approach followed in Ref. E23 is that the trick 
bellow Eq. (3.21) in [22] cannot be used when 11^ is non-zero. 
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a. Scalar modes 


Einstein equations give four equations for the evolution of the scalar modes. The first of 
them comes from <5 Gq = kSTq, and is given by 

* 2 * + 3 H(V + m) - tv [V - 3«] - T_<t, S' - i ik v <r v A - T* [(**)' + 4 H V] 


E 


E' x a T x + 2 E x y VvaVybM^b 

a,b 


2H 11 ” 2 ' 

2«r 

= -a -dp. 


(A39a) 


where, for simplicity, we have introduced the new variable 

A-^ + 4+ (0. 

The second scalar equation can be extracted from 5G ° = k8T® by projecting it along the 
vector k l and is given by 


p(l + w)v + U^v - 


V aVa 


A I a 

(A40) 

The third and fourth equations come from trace and traceless part of 5(S'* = kTT*. They are 

T" + 2'HT / + W + (2W + H 2 )4> - ^ 2 (4> — T) — + ^u 2 (X - 3T) 

3 6 n 11 2 


= a 2 /s; 


^ f -S’a^'tA + 2-E\ y ^va&vbMab J + ^^ ^ ^ + 4'H(J 2 ] 

A \ a,b / a 




(A41a) 


and 

2 


-/C 2 ($ — T) — (Ty 
O 


3 n 




- 2i k y u va $ a - 2A (o') || - — ((cr'),) + 4 y a wa a Yb E x M 


<A 

l a6 


a,6,A 


= a 2 /s; 


+ 57 n nVl + oT7 ES n 


2 T 


3 ^ V a V u ^ 

a A 


2 

- v — y o- TA n TA + 3 y ^n 1 


(A42a) 


respectively. Note that, despite the appearence of i factors, these equations are real. 
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b. Vector modes 


The two equations for the vector modes can be obtained through the combinations 
e a(SG® — k5T®) = 0 and — kST?) = 0. They are given respectively by 


- y<7 va X + j ^ E X Cr vb M 


k 

—2 a 2 n 
k 2 


A 

ab 


b, X 


1 x , T 

p(l + w)v a + ikIi Ya v - -n v a + > M ab Il T \V b - ik—U va 

b, A 


(A43a) 


and 


4>: + + £ i b a TX M x al - A v „A" 


b, A 


+ 


2iT 


kM 

4i 

¥ 


3(J || (o’Ova - 3(T va ((T / ) || +2 Mlb (o- v6 (cr , ) T A ~ <Ga(>')v&) 

b, A 


+ 77^-^ A +-Mife <J vfc( <J T+^A _<T tX^a)] + 7~ ^ EWvbM 


b, A 

~AV) v „ 


A 

ab 


b, A 


2iT 

1VH 


3(7 


(O')va)' - ■y( < 7 , )v“ “ 2 ^K)v<* + Ma,{(r\ h 


b, A 


2i 2 

= —an 
k 




i^r + 2Tn va + 2— a ll n va - —cTyo-cr,! + y] cx vfe n TA .M 


A 

ab 


b, A 


. (A44a) 


c. Tensor modes 

There is only one dynamical equation for the tensor modes. This equation follows from 
the projection (<5G* — rnfTp = 0, which gives 

E a + 2 7-LE' x + k 2 E x — 2E X ^ cr 2 a — 2 A a ct^( 1 _ A ) + 2£ , (i_a)O t+ o tX 


-O'tA 


= a 2 /t 


*U^I+A" 




2X(<j') tX + 2 ifc ^ a v b$ a A 4 ^ - — ((cr') T a)' 


a,6 




2vrJ + n^A — ^ cr N n T A + cr TA n, + y ^ n va <7 v ;,.M 


A 

ab 


a,b 


(A45a) 


4. Perturbed Fluid equations 

The perturbed conservation equation follows from 

{8V fl )TZ + V li (6Tp)=0. (A46) 
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Working in the gauge (A36), the time component (u = 0) of the above expression gives the 
perturbed continuity equation 


Sp + p( 1 + w)W 2 v + 3 U{5p + SP) - (p + P)3V' = 

d j Sir j0 - + HIlijti* - a)5n in lj + a)U u h lj - W j , (A 17a) 


where we remind that h %3 was introduced in Eq. (A19). Likewise, the perturbed Euler 
equation follows from spatial part (u = i ) of Eq. (A46). We hnd 


J^{[p( 1 + w)% i:j + Tl^v 3 } + didP + d l 67r H + (1 + w)pdi$ + II {dfi 

+m{[p(l + w) li3 + II,,; A} - U u djh lj + ^TljdX - l -W k d t h 3k = 0 . (A48a) 


Despite their generality, the above equations are not very useful since they are imple¬ 
mented in real space. In order to obtain their Fourier counterparts we need to project these 
equations along the scalar (k l ) and vector (e l a ) modes. This mode extraction procedure 
is tedious but straightforward, and requires special attention to the use of the evolution 
Eqs. (Al4aa- A14 c|j) of the Fourier wave vectors. 


a. Scalar modes 

Both continuity and Euler equations lead to conservation equations for the scalar modes. 
They are given respectively by 


S' + 3 n[(c 2 s - w)5 + uT] - (1 + w){k 2 v + 3T') - -a ij Ti i3 = - 


P 


P 






nr ij u [<TijJ 


u ij 




+2a 2 K-n ,:? n ? ;, - 8TITV,, + k 2 U^v - i k Y n v0 (u a - $ 0 ) - 6?Dr T + 2k 2 ir s {H + a,,) 


-2iAs Y (Tya^a ~ 2 Y 


^t-Wa 


X] + c7„n TA ) + n A (2 HE X - E' x ) 1 (A49a) 

A A J 


and 


v' + H(1 — 3c 2 s )v + <3> + 


1 + w 


c 2 S + wT — (n 1 — k 2 n b ) 
P 


2i 

a v aVa 


£■ 


a i m i3 (1 + c 


p \ 1 + w 


+ . u ^ p | ifc ( n n v )' + i^ n n $ + + Y n v « {l a) + i^n,, ^ - 1^ - ik Y ^T^ E - 


ik(l + w)p I 

2 Y n v aA« + n®v, 


-JT 


+ Y ( n vA)' + v aV a + 2ika Y Jl va t 


Y n }i a va V a + ?Y a v* V l> M a b n xA ) > , 

a \ b, A ) ) 


(A50a) 



























where we have made use of the equation 


w 


/ 


3n(i + w) + 


n 

p 


w). 


and of the definition [2] 

8 = —, 8P = c 2 s 5p + wpT , 

P 

where c 2 and T are the sound speed and the entropy perturbation, respectively. 


(A51) 


(A52) 


b. Vector Modes 


There is only one conservation equation for the vector modes, which follows from the 


vector projection of Eq. (A50). This equation is given by 


v'a + 'Hil - 3 c 2 s )v a 


k 2 TT^ flb(Tij (l + 


p(l + u;) p \1 + u) J a 2 


] a \\ V a + J2 M afP' l 


A V b = 


b, A 


(1 +U)p 


d v (ikU ya V ^ ^ ^ab^b) ^ ^ ^ab^bc^c ~l~ 4T - /(iA X\.^ a V ^ ^ ^ ab^b) 


b,c 


i/TI va $ + i k I 2 U T xa v bMg b + 2fl va cr || - IL^cr va - %II va 


A, b 


T 

H 


(A53a) 


where we have introduced 


k^ab — ^ij^a^b > ~^ab — 


b ■ 


(A54) 


c. Friedmannian limit 


It is a straightforward exercise to verify that the above equations have a well-defined 
Friedman-Lemaitre limit. Redefining —y P8 ti 1] and ir T —> k 2 n s /3 to compare with 

Ref. [2], we hnd immediately that 


8' + 3n[{c 2 s - to)8 + uT] 
v ' + H( 1 — 3 c 2 )v + <f> 

v'a + W - 3 C 2 )v a 


(1 + u)(k 2 v + W) 

-- — 

1 + U) 1 + u 

u P V 

l + u a ’ 


T-|(2 k\ s ) 


(A55a) 

(A55b) 

(A55c) 


which are the expected equations. 


Appendix B: Perturbed geometric quantities 

This section gathers the expression of the geometrical quantities at order {1,1}, as needed 
for the computation of this article. 
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1. Connections 


Using the commutators of the tetrad field 

7 -gb = ©bf , leab = Vcdl-ab ( B1 ) 

the Ricci rotation coefficients are obtained through 

b Jgbc — Ibd® ) ^ abc i labc R leba I bca ) ■ (B2) 

Up to order {1,1}, the commutators are 

7“oo = 7“oo = 0 ) (B3a) 

7°oi = — 7°io = dft = ~ Pftft , (B3b) 

7 1 = -'ll = -E’j + ftA + Si Bj - a a (l - t) + , (B3 c) 

7 %_ = -2S, ,E j]t + + 2a,,*4 - 2/3’3,'M‘j , (B3d) 

and the Ricci rotation coefficients are thus 

u m = -u m ^ -exjp[-j3\ ij d j ^ , (B4a) 

^ioj — — Ujjo = — <Ty( 1 — $) — Ub + , (B4b) 

^Oij ^ 0 ji 7>?;j , (B4c) 

Wgfc = -cUiAy = 25^8^ + 2—^-2 <%£ fe]i . (B4d) 


2. Riemann and Ricci tensors 


We report the Riemann and Ricci tensor components of the metric (4.1), where the overall 
scale factor a 2 has been removed by a conformal transformation, up to order {1,1}. We first 
give their components in the coordinated basis (with the use of the package xPand [88]). 
and then in the tetrad basis {©}. 

In the coordinate basis, the non-vanishing components are given by 


RoiOj 


Rq ij k 


Rijpq 


lift" + (^)" + raft' + 2*0 - <4(1 - 24/) 
+d,dft - E% + 8 {% B' 3) 


27 njCTkfdq'S' - 2 ai[jdk](p - 2 ( \i\j + dk]^ 

+did\jBk] + 2 d[kE'j^ 


- 47 [i\p<7j} q ft' + 4 (t[Rp + ^r) dj]d q ft 

-dqdjEip - did p E jq + 8 v 8jE iq + 8 q diE jp 


(B5a) 

(B5b) 

(B5c) 
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and 


R 00 = 3$" + 7 vdjdi® 


4cb$ 


Rot = 2di*' - ajdji® + 3$) - -A B t - EJ- 

Rio = 4(1 -2<I>- 2$) - + 3$') + 7 ij [l kq d q d,M> - $"] 

+didj(V - $) + E’r - A Etj - 


n 


+ ^f d k dk^- 2°±d j) d k *. 


Projecting using the tetrad (|4.6|) leads to the components in the tetrad basis 
Ro^ 

+d,d^ - E’r + d (l B’ 3) 


and 


KKOi = 4 $" + ( + 2$') - 4(1 - 2$) 


Rom = 2Siya^Wg^ - 2a i[l dk ] ip - 2 % 4 $' - 2 

+did[jBk] + 2 44 ]* 

Rjjpq 4ft[;[p('Tj,'7]$ T 4 ^4[^[pP(T[/[p/74^ 4] d q j $ 
—dgdjEip - did p E jq + dpdjE iq + d q diE jv 


(Ti [l dk]'&\' 

n 


i?00 = 3$" + 4%$ 


Roi = 2 4$ 7 - 44$ + 3$) - -AS, 


44 vl/ 

R 


% = 4(1 - 2$) - ^($' + 3$') + % 44$ - $"] 
+ 44 ($ - $) + 4 - A4 - 44) 






k 

+ jj-d k d k q< - 2^44$ . 


(B6a) 

(B6b) 


(B6c) 


(B7a) 

(B7b) 

(B7c) 

(B8a) 

(B8b) 

(B8c) 


Appendix C: Glimpse on the full lensing method 

This section details how a tensor field on the sphere X is lensed by a vector field ot, the 
lensing being defined as the result of a parallel transport with respect to this vector field. 

First, for any direction on the sphere, there exists a rotation which connects the azimuthal 
direction with this particular direction n. If this direction has spherical coordinates ( 9 , ip) 
this is simply 


n = R n ■ e z R n = R(tp,9,0) = R z (ip) ■ R y (9) ■ R z (0) (Cl) 
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where R(a,/3,^f) is a general rotation parameterized by Euler angles. Now, if we want to 
define the helicity basis at a given direction n as a result of this rotation applied to the 
hclicity basis at the north pole, we have to face the fact that the helicity basis at the north 
pole is not well defined, since e v is not defined at this point. We choose that at the north 
pole n^ej.) = ^(e x =)=i e y ), since this ensures that the hclicity basis at any point is obtained 
from the one at the north pole through a rotation, that is 

n± ( n ) = Rn ' n± ( e z) • (C2) 

A spin s tensor is defined as X(n) = X s {n)m s {n) = [A" s m s ](n). Its component on the 
polarization basis are simply obtained by projection 

X s (n) = m~ s (n) ■ X(n) = m~ s (e z ) ■ [R-\X(n)} = m~ s (e z ) ■ [R~ l X]{e z ) . (C3) 


This means that instead of projecting a tensor on the polarization basis at a point n we can 
equivalently rotate it so that the point which initially in n becomes located on the azimuthal 
direction. Then we can evaluate its components on the polarization basis at this azimuthal 
direction. The azimuthal direction can thus be used as a common reference for all points 
on the sphere since for each point there is a unique natural rotation to transport from this 
point to the azimuthal direction. 

Let us consider that, due to the lensing vector at, the tensor field we observe in the 
direction n, is now the result of a parallel transport of the underlying tensor by this vector 
field a. Such a parallel transport is equivalent to a rotation around the axis nxa, so that 
the lensed tensor field is related to the unlensed one by 

X(n) = KxocAK”) ■ (C4) 

We use the notation Rv to indicate the rotation defined by the rotation vector V. This 
is the rotation around the axis defined by the vector V with an angle obtained from the 
norm of V. It must not be confused with the previous notation R ni which is the rotation 
that brings the azimuthal direction toward the direction n. As emphasized previously, the 
components of the lensed tensor field, as any tensor field, can be obtained by transportation 
to the azimuthal direction, that is 

X‘(n) = m-(n).X{n) = m-(e z ) ■ ](ej . (C5) 


Using the general property of rotations R n Ry_R n l = Rn n v leads to 


R n R nxa(n) Rn ~ ^e,x[R~ 1 a(n)] ) 


that is to 


D— In-1 _ C> —1 73— 1 

lx n Jl nxa(n) e z x[H~ 1 a(n)] " ’ 


which can be used to recast Eq. (C5) as 


X-(n) = m--(c,).[ir 1 ^ ^(c,). 


(C6) 

(C7) 

(C8) 


This can be understood easily once we extract the helicity components of the lensing vector. 
Indeed, the hclicity basis components of the lensing vector field are obtained just like for 
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any vector field as a ± (n) = n T (n).ct(n ) = e x =F ie y ).[R n 1 ct(n)}. If we define 


Oix = ^(a+ + a_), a y = -J=(a+ - a_), 
= ag cos = ag sin , 


(C9a) 

(C9b) 


where (a x , a y ) are the components of the lensing field once transported to the azimuthal 
direction, and (ag,a,p) their associated polar components, we obtain that 


^'e z x[R n 1 a(n)\ R( a ipi a d, a ip) ■ 


(CIO) 


This means that instead of lensing the tensor held at the point n and subsequently ex¬ 
tracting the component, it is equivalent to transport both the held and the lensing vector 
at the azimuthal direction with R n l and then let the transported lensing vector act on the 
transported tensor held. This procedure is explained graphically in Fig. [ 3 j 



FIG. 3: The azimuthal point Z is at the azimuthal direction e z , while the point P is located in 
R n e z . Because of lensing the signal observed in P is coming from P' whose location is at R nX a(n) • 


j, which corresponds to a covariant transport along a from the point P. As for the point Z', 

e z . This is 


-l 


R r 

it is obtained by applying f?” 1 on P' meaning that it is located at R, 
equivalent to i? e x r R^ l 1 a (n)] ' e z which corresponds to a covariant transport from the point Z, along 


^nxa(n) ' Rn 


the vector R n l (y.{n) 


The lensed field at P (Xp) is obtained by covariantly transporting back along 

— 1 v _ ( r>— 1 

'nx&( 


a the unlensed field at P' (. Xpi ) by application of R ^ , and we get Xp = [ R ‘ , . ■ X 

' ' x ol (ti j y x ol (ti j j ^ 

and it is in general different from the unlensed field at that point Xp. In order to read the 
components, everything is transported back into the azimuthal region by application of R^ 1 - We 


K 1 


R 


X 


For 


get X z > = (R n L ■ X) z „ X z = (Rn 1 ■ X) z , and X z = 
the latter, there is an alternative expressions, which corresponds to covariantly transporting back 


R ?‘ X ) Z = 


-1 

nxcx(n) 


the field X z > along R^a, and this leads to X z = R 




e 2 x[i?„ 1 ct(ri)] 


R" 1 • X 


from which the 


components can be read by projection on the local helicity basis associated with e x and e y . This 


is precisely the meaning of Eq. (C8) which is used to compute the components of the lensed field. 


With this crucial result at hand, we just need to compose the rotations of Eq. (C8) in 


order to obtain the components of the lensed field in terms of the multipole components of 
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(Clla) 


the unlensed tensor field. Using the general transformation law (D16[), we get 


X 8 (n) = Y ln"( e z)D e m " m '[R 1 (a ip ,ag : -a (p )}D £ m , rn [R n 1 ]X^ m 

imm'm" 


Et- 1 )' 


tmm' 


^ ~ n e ro-i 

47T 


D -sm\ R ( a ^ a 6^- a ‘p)\ D m’ml R n}X! rn (Cllb) 




(Cllc) 


tmm' 


where it is understood that the components (ag, a<p) correspond to the lensing vector field 
at the position n considered, and these should be obtained from the definitions (C9). 


From this relation between a tensor acting as a source for an observable and its lensed 
version due to the geodesic structure between the source and the observer, it is possible in 
principle, to obtain the correlations functions (see e.g. Ref. [75] for the case of CMB lensing). 
A simplification can be obtained by expanding the spherical harmonics in a small angle 
approximation. Indeed, very close to the azimuthal direction, the spin-weighted spherical 
harmonics are approximated by 


of + i 

ymp,v) - (- 1 


(C12) 


and given that J m+S {x ) behaves like x m+s when x —> 0 there is a natural way to expand 


Eq. (Cll) in powers of the lensing angle. More details can be found in Ref. ra- 


Appendix D: Mathematical tool-box 


1. From Cartesian to spherical derivatives 


In this section Di refers to the c ovaria nt derivative on the unit sphere in Cartesian coor¬ 
dinates, that is Df 3 defined in Eq. ( |4. 10 ). The key relation to derive all the decompositions 
from Cartesian derivatives di to radial and covariant spherical derivatives is the simple re¬ 
lation 

rdiXj = Sij, (DI) 

which is just the statement that the extrinsic curvature on a unit sphere is equal to the 
metric on this sphere. For a scalar, a projected vector and a projected tensor, we have 


diip 


d i B j 


diEj k 


Di<p 


r 


+ 


D,B 1 ~ ^ x-j ~ 

+ ( Bj) tr Xi ~ yBi 
DiEjk | / 771 \ * Ey 771 fji 


(D2a) 

(D2b) 

(D2c) 


where we use the notation X )T = x l diX for the radial derivative. Note that this radial 
derivative and the covariant derivative on the unit sphere Di commute, as they are just the 
geometric versions of derivatives in spherical coordinates. 
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Iterating these relations we obtain for scalar perturbations 


didjip = 2 x^Dj) + 

V T / ,r 

" ry* £ ry» 

P[x J didj(p\ = Di 0 ) 


DiDjip ~ ~ , r< 'Px 

2 ”1” ,rr%i%j H - &ij 




S?S]d k d p d q¥ > = -S k (jDi) ( ^ ) + Sijdk (^)+d l 


<P,r 


DiDjip 


x k S?S]d k d p d q p = 


DiDjip 


+ S«(^ 


msiw# = l S «jD<) (f) + S tj D k (|f) + 


D k DiDjp 


As for vectors and tensor, the useful relations are 

Ad(i.Ei\ 2 E k j 

x l x J d k diE i j = d k diE r --- - H-— 


-Z rj rj D DiD l B r 

x didjB, = -^— 


r 

2 D l Bi 


+ (_B r ),rr + 2 


2SkiE r ~ ^ 

2 E( k xi) 

2 Br 


( Br ), 


(D3a) 

(D3b) 

(D3c) 

(D3d) 

(D3e) 

(D3f) 

(D3g) 

(D4a) 

(D4b) 


y* ry*£j 

Finally, using the fact that the vector modes are transverse and that the tensor modes are 
transverse and traceless, we get 


Bi 

= Bi + XiB r 

(D5a) 

D l Bi 

= —2 B r — r{B r ) r 

(D5b) 

Eij 

= E^ + 2 EpXj) + E r XiXj 

(D5c) 

S ij Eij 

O 

S- 

1 

(D5d) 

D l Ei 

= —3 E r — r(E r ) }r 

(D5e) 

D l Eij 

= -3 Ej-r{Ej) tr . 

(D5f) 


2. Spin-weighted spherical harmonics 

Spin-weighted spherical harmonics are defined in terms of Wigner D -matrices as [§U] 


[of a. 1 

*&.(«, 0) = y^^(-l) m e is ^i ms (a,/3, 7 ), (D 6 a) 

= (Mb) 

= (-l)”* + ‘y ( 7"(a,/?) (D6 c) 

where a, (3 and 7 are the Euler angles. Wigner D-matrices are in turn defined in terms of 
infinitesimal generators of three dimensional rotations as [89j 

D t mi m a = (im 1 \U{R)\em 2 ), where U(a, /3, 7 ) = e ~ lolJz e~ l/3Jy e~ in/Jz . (D7) 
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In the special case in which the direction is aligned with the z-axis we have 


[of 4 - 1 

= (D8) 

Under the parity transformation 

a —)■ a + 7T, (3 —> n — ( 3 , 7 —>■ 7 + 7r 

the spherical harmonics and Wigner /)-mat rices transform as 

Y lm (a,0) -> (-l)'y fm (o,/?), (D9) 

= (-l)' +m 'By_ m ,(o,/?,T). (DIO) 

In particular, it follows that under parity transformation the Spin-weighted spherical har¬ 
monics behaves as 

YL(a,0) ->■ (Dll) 


3. Rotation of fields on the sphere 

The transformation of the spherical harmonics is given by 

[RY em ]{n) = Yt m (R-\n) = =^lW(n)Dl, m (fl) (D12) 

m' m' 

where the first equality is the definition of the transformation of a function on a sphere 
under a rotation. For a scalar field on the sphere, we can deduce the transformation of its 
multipolar components of its expansion in spherical harmonics, 


A (t'i) ^ ] Tfab m (fi), 

£m 


to be 

[RX](n) =X(R-\n) = Y im ,(n)D e m , m (R)X em 

(D13) 

so that 

£mm' 

[RX] M = Y,&„, m (R)X em . 

(D14) 


m 


The rotation of a tensor field on the sphere is very similar. Once it is broken down into 
symmetric traceless tensors then, by using that such tensors are decomposed as X(n) = 
( n)m ±s {n ) = [X ±s m ±s ](n), it can be expanded in Spin-weighted spherical harmonics 
as _ 

x (») = E WW^'W + x iXMm-(n)} . (D15) 

£m 

Under a rotation, it transforms as (see Appendix A of Ref. [86] ) 

|J?X](n) = R.X(R-\n) = V Y&(n)Di,„(R)X£m±‘(n ), (D16) 

£ m m! 

that is 

{RX\Z = D‘ m , m (R) X%. (D17) 

We remark that it is exactly the same transformation law as for scalar fields because we 
have transformed the full tensor field X ±s (n)m ±s (n) and not just its component X ±s (n) 
considered as a scalar function, for which the transformation law is more complicated [86] . 
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4. Wigner 3j Symbols 


The 3j symbols satisfy the following properties 


M <-2 ^3 

m\ m 2 m^ 


<■2 £3 M 

m 2 m 3 mi 


£3 <-1 t 2 

m 3 mi m 2 


= (-l) 

_ ^_j^h+hj+^2 


h+h+h f ^1 c 3 ^2 
mi m 3 m 2 


tl *3 

—mi — m 2 —m 2 


Moreover, they are identically zero whenever any of the following conditions are violated, 
mi + m 2 + m 3 = 0 , \l t - £ 3 \ < 4 < 4 + 4 , {i, j, k} = {1, 2, 3} . 

They are also orthogonal in the sense that 

4 4 M Mi 4 4 
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and that 
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Since Eq. (D19) holds for any set {mi, m'j, m 2 , m 2 }, two important cases follow from this 
expression. First, consider the case where m 2 = —mi and m 2 = —m\. Then, using the 
selection rule mi + m 2 + m = 0, it follows that 
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mi —mi Oil rn\ —rn[ 0 


= h 


mimj 


(D20) 


Second, note that if we further impose that mi = m\ = 0, then 


Another useful expression is 


E< 2< +!) 


<1 <2 0 
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= 8 , 


Mi 


y/ 24 + 1 


(D22) 


A recurrent expression when dealing with deviations of isotropy is the integral of three 
spherical harmonics, also known as the Gaunt integral, and defined as 
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Note that, due to the symmetries of 3j symbols and the properties of Spin-weighted spherical 


harmonics under complex conjugation, the coefficients s C™^ 2m2 defined in Eq. (5.12) satisfy 
the following properties 


spmim2ra3 _ /_ 1 \<i+£2+4s/^mim2m 3 _ 


—m\ —m 2 —m 3 


(D24) 


From the definitions (D23), ( 5.12[ ) and the closure relation of spherical harmonics one can 
also verify that 

(D25) 


an identity which is needed in order to derive Eq. (5.16) 


Let us also define a useful integral for the gradient expansion approach of lensing by 
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(D26) 


where the polarization basis is voluntarily omitted for a simpler notation. It has the useful 
property inherited from Eq. (D24) 

jm 1 m 2 m 3 _ (]\h+l 2+^3 jm 1 m 2 m 3 fD27i 

t\t 2 t 3 v -F 5 £ 172^3 * 7 

Its expression can be found using 


its 
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= ± ,F hhl , I 11 12 ‘ s I (-1) 
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(D28) 
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where, following Ref. [85], we dehned the symbols 


sFu 3 (.2 = o [^1(^1 + 1) + ^2(^2 + 1) ~ + 1)] 


In particular we have 
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